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Abstract 



We study the properties of the Triple Pomeron Vertex in the perturbative QCD 
using the twist expansion method. Such analysis allows us to find the momenta 
configurations preferred by the vertex. When the momentum transfer is zero, the 
dominant contribution in the limit when N c — > oo comes from anticollinear pole. 
This is in agreement with result obtained without expanding, but by direct averaging 
of the Triple Pomeron Vertex over angles. Resulting theta functions show that the 
anticollinear configuration is optimal for the vertex. In the finite N c case the collinear 
term also contributes. Using the Triple Pomeron Vertex we construct a pomeron loop 
and we also consider four gluon propagation between two Triple Pomeron Vertices. 
We apply the Triple Pomeron Vertex to construct the Hamiltonian from which we 
derive the Balitsky- Kovchegov equation for an unintegrated gluon density. 
In order to apply this equation to phenomenology, we apply the Kwiecihski-Martin- 
Stasto model for higher order corrections to a linear part of the Balitsky-Kovchegov 
equation. We introduce the definition of the saturation scale which reflects properties 
of this equation. Finally, we use it for computation of observables, such as the F 2 
structure function and diffractive Higgs boson production cross section. The impact 
of screening corrections on F 2 is negligible, but those effects turn out to be significant 
for diffractive Higgs boson production at LHC. 



Zusammenfassung 



Wir untersuchen die Eigenschaften der Triple-Pomeron- Vertex in der 
storungstheoretischen QCD mittels einer Twistentwicklung. Eine solche Analyse 
ermoglicht die Bestimmung der von der Vertex bevorzugten Impulskonfiguration. 
Es ergibt sich, dass fur N c — > oo und verschwindenden Impulstransfer der do- 
minierende Beitrag von dem antikollinearen Pol stammt. Dies stimmt mit dem 
Ergebnis iiberein, das man erhalt, indem man auf eine Twistentwicklung verzichtet, 
jedoch iiber die Winkel der Vertex mittelt. Die resultierende Theta-Funktion 
zeigt dann, dass die antikollineare Konfiguration bevorzugt ist. Dies andert 
sich fur endliches N c , da in diesem Fall auch die kollinearen Terme beitragen. 
Mittels der Vertex konnen wir einen Pomeron-Loop konstruieren. Dariiber hinaus 
betrachten wir den Fall in dem vier Gluonen zwischen den Vertices propagieren. 
Wir benutzen die Triple-Pomeron- Vertex um die Balitsky-Kovchegov-Gleichung fur 
die unintegrierte Gluonendichte abzuleiten. 

Um diese Gleichung in der Phanomenolgie anwenden zu konnen, benutzen wir 
das KMS-Modell um Korrekturen hoherer Ordnung in den linearen Teil der 
Balitsky-Kovchegov-Gleichung einzuschlieBen. Wir fiihren eine fur diese Gleichung 
natiirliche Definition einer Saturierungsskala ein. Wir benutzen diese schliefilich, 
um Observable wie die Strukturfunktion F 2 und den Wirkungsquerschnitt fur 
diffraktive Higgsbosonproduktion zu berechnen. Es stellt sich heraus, dass die 
Bedeutung der Abschirmeffekte fiir F 2 vernachlssigbar ist, wahrend sie fur die 
diffraktive Higgsproduktion am LHC bedeutsame Beitrage liefern. 
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Chapter 1 
Introduction 



The strong interactions in a field theory language are described by Quantum Chro- 
modynamics (QCD). Quark and gluon fields are its degrees of freedom. Its structure 
and its predictions for phenomenology can be investigated in a perturbative or non- 
perturbative way, depending on the question we ask. In this dissertation, we are 
going to address the questions which are tractable within the perturbative approach, 
but which also shed light in the regime when nonperturbative physics starts to be- 
come important. In particular, we are interested in the high energy limit, or Regge 
limit of QCD, which is defined as the limit where the energy of colliding particles 
is much higher then other scales. To be explicit let us consider two virtual photons 
scattering. The leading order diagram that contributes is just two gluon (in singlet 
state) exchange in the t-channel. To have a realistic description one should, how- 
ever, calculate higher order corrections in order to check the validity of the Born 
term (just two gluons). The higher order corrections, which can be separated into 
a real emission part (gluon in the s-channel) and a virtual emission part (gluon in 
the t-channel) , turn out to be proportional to the logarithm of the energy multiplied 
by a strong coupling constant. Since the energy can be large the smallness of the 
coupling constant may be compensated and the whole series should be resumed. 
The summation of the large logarithms of the energy accompanied by the strong 
coupling constant is performed by the Bethe-Salpeter like equation i.e. BFKL [1] 
(Balitsky, Fadin, Kuraev, Lipatov) integro-differential equation. This equation is the 
basic mathematical tool which enables us to investigate the Regge physics within 
the QCD and it is believed that one should go along the line of BFKL equation to 
understand its high energy behavior. In the color octet channel, BFKL describes 
the propagation of a reggeized gluon. It is the gluon which represents collective 
excitation of the gluonic field and which turns out to be the basic degree of freedom 
in the high energy limit of QCD. On the other hand, in the color singlet channel, 
the BFKL equation defines so called BFKL pomeron which is a bound state of two 
reggeized gluons, and which posses vacuum quantum numbers. Studies show [2] 
that LO BFKL is not enough and higher order contribution is necessary in order 
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to make BFKL applicable to phenomenology. What is even more important is the 
fact that at very high energies BFKL amplitude violates the unitarity bound [3] 
which states that the amplitude of hadronic process should not grow faster than the 
energy logarithm. One cannot saturate this bound using the linear equation (even 
including subleading corrections). The total cross section of a scattering process, as 
obtained from the BFKL, grows like a power of the center of mass energy squared s. 
The corresponding exponent equals at leading order ojbfkl — 0.5. The violation of 
unitarity is related to the fact that at very high energies we apply the BFKL equa- 
tion to physical system which is very dense (for example Deep Inelastic Scattering 
at Bjorken x = 1CT 6 ) and where the probability of recombination, i.e. probability 
for fusion of gluons, is high and this is not included in the BFKL equation. Such 
processes reduce the growth of the partonic density and lead to the effect called 
perturbative saturation. 

There are two main streams of activity to construct unitary theory of high energy 
scattering. One of them is persuaded in momentum space where reggeized gluons are 
its basic degrees of freedom. Here one can distinguish two approaches: one based 
on Lipatov's effective action [4], and the other based on investigations of Bartels 
[5]. The latter approach is of particular interest for us. Bartels [5] generalized the 
BFKL equation to situation where more then two reggeized gluons (in t-channel) 
are exchanged between colliding objects [5, 6, 7]. The physical process that one con- 
siders here is again scattering of two virtual photons. The integral equations which 
generalize the BFKL are written for an amplitude which is a coupling of gluons to 
virtual photon via the quark loop. The solutions of the considered equations are 
used to construct imaginary part of amplitude for photon-photon scattering. This 
approach leads to the emergence of the field theory like structure where transitions 
from n to m reggeized gluons are possible (in contrast to BFKL Green's function 
which is just the propagation of two gluons or Bartels, Kwiecihski, Praszalowicz 
(BKP) [8, 9] which describes propagation of constant number of n gluons). The 
integral equations for n gluon amplitudes form an open set of integral equations. 
The equation for the three gluon amplitude can be solved. The equation for the 
four gluon amplitude can be partially solved. It can be rewritten as an equation 
with reducible and irreducible parts. The basic ingredient of the irreducible part 
is a two-to-four reggeized gluons transition vertex. This vertex plays a central role 
in the physics of saturation because it allows gluons to recombine and therefore to 
tame the strong power like growth of gluon density and therefore the cross section. 
There is a hope that the approach outlined above can be mapped to the conformal 
field theory because of the fact that Mobius symmetry is encoded in the BFKL 
equation and also in the two-to-four transition vertex [10]. If it would be possible 
one could use the methods of conformal field theory to study higher order vertices 
and perhaps to solve this theory. 

Another approach to high density systems is based on the dipole degrees of freedom 
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of the high energy limit of QCD. The basic equation here is again BFKL, but for- 
mulated in the coordinate space [11]. That equation is written for an amplitude for 
scattering of an color dipole off a target (which may be another dipole). The evolu- 
tion of a dipole is defined by summation over a cascade of dipoles scattering off the 
target. There is, however, assumption that every child dipole scatters independently 
of others. When the density of dipoles is high one should allow for simultaneous 
scattering of child dipoles what leads to set of coupled integro-differential equa- 
tions known as Balitsky hierarchy [13]. The equation for scattering of a one dipole 
contains contribution from two dipoles. In turn, the equation for two dipoles, con- 
tains contribution from three dipoles and so on. That chain can be broken under 
assumption of a large target and large N c . In particular, if we are interested in 
the first equations from that hierarchy, we can assume that the amplitude for two 
dipoles scattering is proportional to the product of one dipole scattering. Under this 
assumption nonlinear integro-differential evolution equation known as the Balitsky- 
Kovchegov equation [13, 14] can be obtained. The linear part, where creation of 
gluons is described by the BFKL term, is accompanied by the nonlinear term which 
is given by the Triple Pomeron Vertex (the same as in the t-channel approach) [6] 
and which allows for recombination of gluons. Balitsky-Kovchegov (BK) equation 
leads to the unitarisation of the dipole amplitude at fixed impact parameter. It is 
the best and relatively simple equation which can be applied to phenomenology. 
This equation is, however, theoretically satisfying. The proper evolution equation 
which could describe saturation physics should incorporate pomeron loops. Pomeron 
loops, or, quantum fluctuations of, pomeron are required to ensure projectile tar- 
get symmetry which is violated by BK. They are expected to be important at low 
parton densities and may be important for understanding of double diffractive pro- 
cesses which probe their structure. Another limitation of BK comes from the fact 
that the BFKL kernel and the Triple Pomeron Vertex are considered at leading or- 
der in \nl/x. From studies of Kwiecihski, Martin and Stasto [2] it is known that 
subleading corrections (which are modeled in consistency with exact NLOlnl/x [15] 
calculation) to linear part reduce the value of the BFKL intercept and lower the 
normalization of its solution. Therefore, in order to have equation applicable in 
phenomenology, one should incorporate those corrections at least in linear part of 
BK. The NLO lnl/x corrections to the Triple Pomeron Vertex are still not known. 
It is not clear how to model those corrections in a consistent way. 
Some of the problems of high energy QCD listed above have been already attacked. 
There are promising approaches [16, 17] to include pomeron loops in evolution equa- 
tion. There NLOlnl/x corrections to BFKL are also known[15] so there is hope that 
NLOlnl/x corrections to Triple Pomeron Vertex can be calculated in a similar man- 
ner. What is missing and what we find important from the theoretical point of 
view is the better understanding of properties of the Triple Pomeron Vertex which 
is the basic element of already mentioned pomeron loop and of the BK equation. 
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By knowing the high energy QCD problems at high energy we are going to improve 
the current state of understanding. 

Understanding of the properties of the Triple Pomeron Vertex is the main subject 
of this thesis. We focus in particular on dependence on momentum variables in case 
of strong ordering which is crucial for understanding of the region of phase space 
where the pomeron loops become important. This vertex also appears in perturba- 
tive approach to diffractive processes. Diffraction is still not really understood and 
we hope that our analysis will be useful in understanding it. Another motivation 
to study momentum dependence of the TPV is to clarify relations between different 
nonlinear evolution equations that appear in the literature (GLR, GLR-MQ, BK). 
To study properties of TPV we apply the method of twist expansion which allows 
us to perform the analysis via expanding considered amplitude in terms of a largest 
momentum. We also find it important to derive analogous equations to Balitsky 
equations directly in momentum space. Here the interpretation is much easier and 
follows directly from Feynman diagrams. This is somewhat obscure in the original 
formulation. Later on, assuming mean field approximation, we derive the BK equa- 
tion for the unintegrated gluon density which is in agreement with the one obtained 
from transformation from coordinate space to momentum space. 
As we have already mentioned the BK equation misses important subleading correc- 
tions and therefore can not be directly applied for phenomenology. In order to use 
that equation to describe the data we propose to model the subleading corrections 
in In 1/x following [2]. This approach turned out to be very successful in describing 
F 2 HERA data and is consistent with NLOlnl/x corrections to the BFKL kernel. 
The next important point is the fact that, due to nonlinearity, the gluon density in 
saturation region strongly depends on the profile function of the proton. In order to 
have realistic description of F2, we study the properties of the gluon density in two 
important cases: the cylinder like profile function and the Gaussian-like. We argue 
that the later is more accurate. Finally, equipped with the fitted nonlinear equation 
for the gluon density, we ask ourselves whether saturation effects, which result in 
hard rescattering corrections, may have any impact on diffractive Higgs production 
at LHC. In existing calculations such effects were not considered. In order to es- 
timate such corrections we modify the standard approach introducing rescattering 
corrections via the Triple Pomeron Vertex and the gluon density provided by the 
BK equation. Our estimation is rough, but can be considered as a first step towards 
more accurate description. 
The thesis is organized as follows: 

• chapter two - in this overview chapter we discuss linear evolution equations of 
QCD and give arguments why saturation effects are needed. 

• chapter three - we construct amplitude for elastic scattering of hadronic object 
which exchange pomeron loop and also BKP state. We perform the collinear 
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analysis of the Triple Pomeron Vertex. 

• chapter four - is devoted to derivation of the BK equation for the unintegrated 
gluon density from an effective Hamiltonian consisting of the BFKL kernel and 
the Triple Pomeron Vertex. We also discuss the relation of the BK equation 
to other nonlinear equations. 

• in chapter five we apply the BK equation with subleading corrections in In 1/x 
to describe the F 2 structure function and we calculate a saturation line. We 
also calculate hard rescattering corrections to diffractive Higgs production at 
LHC energies. 

We end this thesis with summary and conclusions. 



7 



Chapter 2 

Evolution equations of QCD 



2.1 Linear evolution equations of QCD 



In this overview chapter we are going to set up conventions that will be used later 
on in this thesis. We start with the Deep Inelastic Scattering of an electron off 
a proton which is a basic scattering process in QCD. Then we introduce DGLAP 
evolution equation which is a realization of the renormalization group in QCD. In 
the next step we introduce the BFKL evolution equation which is the first step in 
perturbative QCD towards understanding of the physics of dense partonic systems. 
Finally we give arguments why rescattering corrections, where the Triple Pomeron 
Vertex plays a central role, are needed, and we give examples of formalisms within 
which those corrections are discussed. 

2.1.1 DIS kinematics and variables 

The high energy limit of QCD can be explored thanks to the asymptotic freedom, 
in a perturbative way. To introduce perturbative parameters and quantities with 
which we will work let us consider the scattering of an electron off a proton at high 
momentum transfer 



where X is an undetected hadronic system. The quantities that are measured are 
the energy and the scattering angle of the outgoing electron. This process can be 
described using standard variables: 



e(p) + N(P)^e'(p')+X(P x ), 



(2.1) 



q 2 = (p-p') 2 = -Q 2 <0, 



(2.2) 



the four momentum squared of the exchanged photon, 



Q 2 Q 2 



(2.3) 



x = 



2p ■ q Q 2 + W 2 - m 



2 ' 

'N 
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Figure 2.1: Deep inelastic scattering process. 



the Bjorken scaling variable and 



W 2 + Q 2 - m\ 



N 



(2.4) 



y = 



the inelasticity. 

In the formulas above W 2 = (q + P) 2 is the energy squared of the photon- nucleon 
system and s = (p + P) 2 is the total energy squared of the electron-nucleon system. 
The nucleon mass is denoted by m^r. When Q 2 involved in the scattering considered 
above is much larger than m 2 N we speak about Deep Inelastic Scattering. 

In perturbative QCD (pQCD) the cross section for process depicted in Fig. 2.1 
is expressed via kinematical variables and dimensionless quantities that contain the 
information about the structure of the nucleon and which are called structure func- 



Airn 1 1 1 

- 2 = -^f foWiO^ Q 2 ) - F*(x, Q 2 )) + «[! + (!- v) 2 ]P2(x, Q 2 )] . (2.5) 



In pQCD structure functions can be calculated in two approaches: 

• collinear factorization where the basic evolution equation in the momentum 
variable Q 2 is the DGLAP equation 

• high energy factorization where the basic evolution equation in the energy 
variable s is the BFKL equation 

In the two following sections we are going to give an overview of those approaches. 



In the collinear factorization approach the structure function is written as a convo- 
lution in longitudinal momentum fraction: 



tionsFx^Q 2 ), F 2 (x,Q 2 ): 




2.1.2 DGLAP evolution 




(2.6) 
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Figure 2.2: Ladder diagrams contributing to DGLAP evolution. 



where Ci(y) (which is calculable in a perturbative way) is a coefficient function 
which represents the hard subprocess cross section for an incoming parton. At 



leading order (LO) in a s , C q = e 2 5(l 



z), C- q 



e 2 6(l - z), C g = 0. The other 



quantity xfi(x,Q 2 ) is the parton density of type i and represents the probability 
to pick a parton of type % with longitudinal momentum fraction x and virtuality 
smaller than Q 2 out of the hadron. To relate parton densities at different scales one 
uses the DGLAP evolution equation which is a realization of the renormalization 
group in QCD. To introduce that framework let us consider scattering of a highly 
virtual photon with a virtuality Q 2 on a proton which is characterized by a scale 
Qq<:Q 2 . According to the naive parton model when the structure of the proton is 
resolved, the probe interacts with point like partons inside the proton. The parton 
density turns out to be independent on the momentum scale and depends only on 
x. This is known as Bjorken scaling. However, if we introduce QCD corrections, 
which allow for parton interaction, the situation changes. If we change (increase) 
the virtuality of the probe we realize that the parton itself is surrounded by a cloud 
of partons, so consequently the parton density grows effectively towards larger Q 2 . 
The dependence of the parton density on the momentum scale can be written as 
integro-differential evolution equation: 



dxf t (x,Q 2 ) _ a s (Q 2 
<91nQ 2 /Q2 " 



2tt 



x 



dy X -f PrMj) 



(2.7) 



The Pij terms are the so called splitting functions which describe the probability of 
finding a parton with momentum fraction x within a parton of momentum fraction 
y. The scale Ql acts as a cutoff, above which we may use perturbative methods. 
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Equation (2.7) can be interpreted as summation of ladder diagrams (in a physical 
gauge) (2.2) of whose rungs are strongly ordered in the transverse momentum. The 
natural evolution parameter in (2.7) is the resolution lnQ 2 / m< 3o (we shall put 
Ql = lGeV 2 from now on). The full set of DGLAP equations describing evolution 
of quarks and gluons takes the following form: 



dxg(x,Q 2 ) a s (Q 2 



d\nQ 2 2tt 

dZ(x,Q 2 ) _ a s (Q 2 ) 
d\nQ 2 2tt 



dy 



dy 



p qq ( y )-x - q 2 +p q9 (y,Q )-9 ( -,Q 
y \y J y \y 



Qns(x,Q 2 ) a s {Q 2 ) 



^£ fx r 
I dyP qq (y)-q NS l-,Q' 

Jx y \y 



(2.8) 



(2.9) 



(2.10) 



<91nQ 2 2tt 

where S(x, Q 2 )=^2i[qi(x, Q 2 ) + <?«( a; > Q 2 )} is the singlet quark distribution and q^s = 
q(x, Q 2 ) — q(x, Q 2 ) is the nonsinglet quark distribution. The splitting functions read: 



P qq {x) = C F 



T R [x 2 



P gq (x) = C F 



1 + x 2 
1-x 2 ). 

(l-x)% T R 

l + (l-:r) 2 

x 



1 



(2.11) 

(2.12) 
(2.13) 



where 



P 99 {x) = 2C A 



x 



1 - X /-, x 

+ + x(l-x) 



(1 — x) + X 
+ 5(l-x)^ CA - 6 An ^ 



fix) 



o ^ (1 - x) + 







X 



(2.14) 
(2.15) 

(2.16) 



and 1/(1 — x)+ — 1/(1 — x) for x < 1. 

At very small x, which is our main interest, the dominant role is played by the gluons. 
This can be observed analyzing the low x limit of splitting functions. The splitting 
function P gg gives the most dominant contribution and may be approximated by 
Pgg^ 2 ^- Neglecting the quark contribution which is subleading at low x we obtain 
from (2.8): 



dxg(x,Q 2 ) a s (Q 2 



^£ f X r 

/ dyP gg (y)-xg I -,Q S 
j x y \y 



i)\n(P 2tt ; ' " \ ■ I ; --' 7! 

This equation can be solved analytically in DLLA (Double Leading Log Approxi- 
mation). Taking the low x limit of P gg and fixed coupling constant (calculation with 



11 



running coupling constant does not change the result very much) and using a saddle 
point approximation we obtain: 

xg(x, Q 2 ) = xg(x, Cfi) expJ^-\n %. (2.18) 

Here we reintroduced Ql for a clearer presentation. The gluon density (2.18) grows 
fast with energy and when used for calculation of cross sections may violate the 
unitarity constraints at small x (more detailed discussion about unitarity will be 
presented later). 



2.1.3 BFKL approach 

An alternative approach to the evolution in the momentum variable (more suitable 
for high energy processes) is to consider fixed momentum variable and to perform 
evolution in the energy variable. Such an evolution is achieved by the BFKL equation 
which can be formulated as a Bethe-Salpeter equation. The perturbative terms 
that are summed up are logarithms of energy accompanied by the strong coupling 
constant (a s In s/s ) n , s ~M hadron (lns/s = lnl/x). In order to justify perturbative 
approach the coupling constant has to be small. The dedicated process where one 
could expect the BFKL resummation to appear is the scattering of two virtual 
photons at high center of mass energy y/s. The BFKL integral equation reads: 

(u - u(k) - u(k + q))d 2)K * <} (k, k', q) = 0( 2 )°{*.«0( k , k', q) 

' n ' tf&- W (l, q - 1; k, q - k)^ 2){b "<>(l, k', q) (2.19) 



/i 



(27r) 3 k 2 (k-q) 2 

where k',k,q, 1 are two dimensional vectors. They represent the transverse compo- 
nents of momenta of the gluons that live in a plane transverse to the collision axis. 
The transverse momenta are not ordered (in contrast to DGLAP), but are of the 
same order around sq. The longitudinal components are However, the longitudinal 
components are ordered what leads to ordering in rapidity and gives rise to BFKL 
summation. 
In (2.19) 

is the BFKL Green's function and p a ^ a 2, a 'i a ' 2 i s the color projector. We are here in 
particular interested in a color singlet channel such that the color projector is given 
by: 

Xaia2Xa' 1 a' 2 

pa 1 a 2 ,a' 1 a 2 _ u In cy\\ 

iV 2 — 1 y ' ' 

In (2.19) the index i runs from 1 to 2, and u is a variable conjugated to the rapidity 
Y via the Mellin transform: 

(Y; k, k', q) = J ^-e Y ^ (k> ^ q) (2 . 22) 
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The inhomogeneous term is given by: 



£(2)0K,aU( k)k ' i)r ) = ( 2?r )3 (5(2) ( kl k i) pKX> (2.23) 

k^ki-r) 2 



and 

2 k 2 (l-q) 2 Z 2 (k-q) 



K {a}->{b} _ 2f f 
Iy 2^2 ~ y J haicj ca 2 b 2 



q 



(2.24) 



(k - l) 2 (k - l) 2 

is the real emission part of the BFKL kernel which mediates transition from two-to- 
two reggeized gluons. The term: 

-" k » = -^/ ^ + »-DMk-i)' (2 ' 25) 

is the virtual part of the BFKL kernel which cancels the singularity of the real part 
arising at k = 1. The cross section for scattering of hadronic objects within the 
BFKL formalism can be written in a factorized form: 

°{Y,Q 2 ) = 2ir J (2^3 l^ foi. Q)G matM (Y, kx, k 2 , q = 0)# } (k 2 , Q), 

(2.26) 

where the functions 1 tia2 =5 aia2 1 , 02 1<l2 =<5 aia2 02 are the impact factors and represent 
coupling of external particles to the Green's function. Indices a iy bi refer to the 
color degrees of freedom. If we are in particular interested in DIS processes we 
define the unintegrated gluon density considering the scattering of a virtual photon 
</>t,l (ki, Q)(transversally or longitudinally polarized) on a proton and setting in 
(2.26) as the 2 proton impact factor and rewriting it as 

<Jt,l{Y,Q) = J ^0T,L(k 1 ,Q)/(F,k 2 ), (2.27) 

where 

k2) = / 0f k4<f)p{k ' )g{Y > k ' k ' } (2 ' 28) 

is the unintegrated gluon density. The color factor Cf, resulting from contracting 
color indices, has been absorbed in the definition of the proton impact factor. From 
the point of view of QCD the proton impact factor is of nonperturbative origin and 
can only be modeled. The virtual photon impact factor vanishes when momentum 
variable equals zero, indicating infrared finites. 
The usual gluon density is then given by: 

r Q2 r/k 2 

xg(x,Q 2 ) = J -^f(x,k 2 ) (2.29) 
and the structure functions by: 

f)2 r J2l, 

F 2 (x, Q 2 ) = / ^<Mk, Q)f(x, k 2 ), (2.30) 
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Figure 2.3: High energy factorization applied to DIS. 



where 0(k, Q) = (pxik, Q) + 0i(k, Q) and we found convenient to use x instead of 
Y— log(l/a;) in order to use the standard DIS DGLAP notation. Similarly, as in the 
previous section we may ask what is a prediction for the gluon distribution in the 
unintegrated case. To answer that question we write the BFKL equation directly as 
an evolution equation in In 1/x and since we are interested in the leading behavior, 
we consider scattering in the forward direction: 



agjx, k, k' 

d In 1 jx 



N c a s 2 



7T 



dn 2 



K, 



g( x , k 2 , k' 2 ) - g( x , k 2 , k") g( x , k\ k 



I2\ 



+ 



^4k 4 + k 4 



(2.31) 

and using relation (2.28) we write the BFKL equation for the unintegrated gluon 
density as: 



df(x,k 2 ) 
d In 1/x 



7T 



dn 2 



f(x,K 2 )-f(x,k 2 



f(x,k 2 



'4k 4 + k 4 



(2.32) 



Such an equation can be solved taking the Mellin transform with respect to the 
momentum variable. The procedure leads to diagonalization of the kernel. The next 
step is to solve the differential equation in rapidity In 1/x and to invert the Mellin 
transform. To invert the Mellin transform the leading eigenvalue of the kernel is 
taken and integration is performed in the saddle point approximation yielding: 



f(x,k 2 



x 



-X r 



ln(xo/x) 



-.1/2 



exp 



2\ 1 



C4 \n(x /x) 



(2.33) 



with A = ^^41n2 w 0.5 for a s = 0.2. Xq is the value of x at which the evolution 
starts and C = 3a s 14£(3)/7r. We see that the dominant contribution is given by 
the first term and the predicted growth of the gluon density is very large. When we 
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apply that result to estimate photon-proton cross section (2.27), the x dependence 
of it is: 



A / \ A 
X \ IS 



That result is in contradiction with Froissart-Martin theorem stating that hadronic 
cross section should behave as: 

ctot < const In 2 s (2.35) 

The authors of that theorem required that the S matrix, containing the information 
on the scattering process, has to be unitary and that the range of interactions 
is finite. The other terms in (2.33) reflect the diffusion properties of the BFKL 
equation. One sees that f(x, k 2 ) is a Gaussian distribution in ln(k /k ), with a 
width growing as ^/\n(x /x) when x — > 0. The diffusion property of the BFKL 
causes problems due to the possible diffusion into the nonperturbative regime, where 
the perturbative approach is not reliable. To cure those problems and to describe 
physics at low values of x, a generalization of the BFKL equation is needed. The 
goal is to have an evolution in rapidity, so that the essential BFKL approach is 
preserved. The generalization can be understood in many ways. In this thesis we 
are going to consider the so called rescattering corrections to the BFKL equation 
which introduces nonlinear effects into the evolution process. This is the simplest 
modification one can do. The outline of the approaches which are in line of the one 
presented in this thesis will be sketched in the next section. 

It is of particular interest to note here that (2.32) coincides in the collinear limit 
(defined as k 2 >K 2 ) with the DGLAP equation for xg(x,Q 2 ), when the dominant 
contribution to P gg at low x is taken. To be explicit we observe that taking collinear 
limit in BFKL we get: 



df(x,k 2 ) N c a s f k d* 2 



f(x,K 2 ) (2.36) 



dlnl/x it J k 2 
Integrating over ln/c 2 (we set k 2 = lGeV 2 ) and using (2.29) we obtain 

dxg(x, k 2 ) N c a s 



<91nl/x<91nk 2 7i 



-xg(x } k 2 ) (2.37) 



what coincides with (2.17) when P gg (the dominant contribution to it) is inserted 
and differentiated with respect to Inl/x. 



2.2 Rescattering corrections 

The approaches described above which are based on linear evolution equations, 
do not take into account so called rescattering corrections. In case of DGLAP 
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one usually says that it takes into account so called leading twist contribution to 
the structure function. This means that in the DIS case the contribution to F 2 is 
assumed to come from two gluons in the t-channel. In case of the BFKL equation 
one should keep in mind that the gluon is in fact reggeized what means that it gets 
contribution from all twists. However, it does not include rescattering corrections 
which manifest themselves in high energy limit of QCD as a propagation of more 
than two reggeized gluons in the t-channel. The physical picture behind the need 
for including those effects follows when one considers parton system at small x. 
At small x the system of partons becomes so dense that partons start to overlap 
and screen each other. Thus we have to account for recombination. Those effects 
are often called rescattering corrections. The way to describe them in momentum 
space in terms of Feynmann diagrams is to consider more than two gluons in the 
t-channel which goes beyond the leading order approximation. In particular, within 
the BFKL framework, the program to include systematic higher twist effects has 
been formulated by Bartels and the basic elements of hopefully underlying field 
theory (conformal), have been constructed by Bartels [5] and Wiisthoff [6], Bartels 
and Ewerz [7]. This approach is usually called the Extended Generalized Leading 
Log Approximation (EGLLA). 

Another approach within pQCD to the high energy evolution and in particular to 
the unitarisation problem is addressed in configuration space. It is also possible 
to formulate in coordinate space the high energy scattering in terms of dipoles 
degrees of freedom. The basic equation, as far as one will not address the problem 
of unitarisation, is the dipole version of the BFKL equation. However, when one 
allows for evolution to the higher energies, one has to include additional processes 
which are modeled by the BK equation. Below we give more detailed outline of 
those two streams of activity. 

2.2.1 EGLLA 

The main objects of considerations in this approach are n (reggeized) gluons am- 
plitudes for which coupled integral equations have been formulated [7]. The gluons 
are on the projectile side coupled to the photon impact factor and propagate in 
the t-channel. The coupling from the target side is not fixed. We are in particular 
interested in the Triple Pomeron Vertex which emerges at the level of four gluons. 
The relevant equations up two four gluon amplitude are given by: 

2 

(uj-J2 "WW" = £> 2 T + *fir {o} ® D h J b \ (2.38) 

3 

(to - u(k t ))D a 3 ia2a:i = D™ 3 + K^ {a} <g> D b 2 lb2 + K 2 b l7 {a} ® D b ^ b \ 
i=i 

(2.39) 
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(u -J2 iU (^)) D T a2a3a4 = D a ^o 2 + K%£ {a} <g> D b 2 lb2 + J2 K 2 b l7 {a} ® D b 3 lb2h 
i=i 

+J2 K ^ {a} ® Db i b2Hb4 i 

(2.40) 

The first of those equations is identical with the BFKL equation. The inhomoge- 
neous term D 2 }q 2 = <f) aia2 is given by the quark loop. The superscripts refer to a 
color degrees of freedom of reggeized gluons. The convolution symbol <E> stands for 
an integral over the loop momentum and propagators p- for each of two gluons 
entering the kernel. 

The novel term appearing in (2.39) is the kernel K 2 ^^ a ^ which mediates tran- 
sition from two to three gluons. 

It's structure is similar to the BFKL kernel and reads: 



7^{6}->{a} _ 3 f r f 

xl 2^3 iJ Jh\a\cj ca^dj dao,b2 



:kl + k 2 + k 3 r ^^^^ ^^^^ 



(k 3 - q 2 ) 2 ( kl - Ql )2 

4 



+ 



q 2 q|k 2 



(ki- qi ) 2 (k 3 -q 2 ) 2 . 

(2.41) 



The kernel K^}^^ which mediates transition from 2 to 4 reggeized gluons reads: 

q2( kl + k 2 + k 3 ) 2 



7^{&}-»{a} _ 4 r f f f 
xl 2^4 y JbiaicJ ca^dj da-jej ea.462 



(ki + k 2 + k 3 + k 4 ) 2 - 



(k 4 - q 2 



,2 



q 2 (k 2 + k 3 + k 4 ) 2 + q 2 q 2 (k 2 + k 2 ) 



(k!- qi )2 ( kl _ qi )2 (k4 _ q2) ; 

(2.42) 

The equation (2.39) can be solved with an initial condition provided by the quark 
loop with attached three gluons: 

D a ^ 2WA = \ fa ia2 a 3 [£> 2; o(ki + k 2 , k 3 ) - L> 2; o(ki + k 3 , k 2 ) + ^(kx, k 2 + k 3 )] (2.43) 

The expression above follows from reggeization property of the gluon at high ener- 
gies. It states that coupling of three gluons to the quark loop can be expressed as 
combination of couplings of two gluons to the quark loop. The solution turns out 
to be expressible in the same way in terms of two gluon amplitude: 

D a ia2 a 3 (ki) k2 ks) = ^ gfam [£> 2 ( 12> 3) - D 2 (13, 2) + D 2 (l, 23)] (2.44) 

where we have applied the shorthand notation D 2 (ki + k 2 , k 3 ) = D 2 (12, 3) etc. The 
structure of the solution expresses the fact that a compound state of three reggeized 
gluons does not exist (at least as one considers their coupling to the quark loop). 
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The equation for four gluons cannot be completely solved. It can, however, be split 
into a reggeizing part D± and an irreducible one with respect to reggeization D\ 

jja\a2aian jjR\aia2aza^ _|_ ^7501020304 ^ 

where 

D^ 1020304 (ki,k 2 ,k3,k4) = -g 2 d aia2a3a4 [D 2 (123,4) + £> 2 (1,234) - L> 2 (14,23)] 
-g 2 d ama * [D 2 {IU, 2) + ^2(124, 3) - L> 2 (12, 34) - £> 2 (13, 24)], 

anrl ,701020304 1 X X 1 1/J ,7 f f N 

cuiu U 2 JV C 0i02 u ai02 ~r 4 l u ai a^k^ka^a^ J a\a2kjkao,at t ) 

For £)^ aia2a 3 a4 a new integral equation can be derived: 



[uo- 



-^^(A: l )) J D 4 /aia2asa4 (k 1 ,k 2 ,k3,k 4 ) = v 2 a ^ ;aia2a3a4 ®^^+^xg^ {a} ®L>f 1626364 



(2-46) 

In the above expression y^^°^ i s an effective vertex which mediates transition 
from two to four reggeized gluons. It represents a compound state of four gluons 
and cannot be reduced to a sum of two gluons states. It, however, can be expressed 
in terms of the so called G function which has the momentum structure as the 
kernel K^}^^ plus trajectory functions. The explicit formula for that vertex will 
be presented in the next chapter. This vertex has the property of being gauge 
invariant i.e. it vanishes if any of the momentum arguments is zero. It is also 
symmetric with respect to the simultaneous interchange of color and momentum 
indices. A further property of the vertex (as in the case of the BFKL) is its invariance 
under Mobius transformations. The two-to-four transition vertex plays a crucial role 
in construction of theoretical framework for saturation physics because it allows to 
include rescattering effects. 



2.2.2 Dipole approach 

The alternative approach to unitarisation of high energy QCD is persuaded in so 
called dipole picture [11, 38]. In this framework one considers decomposition of the 
photon wave function into its Fock components. In the lowest order, this is the bare 
quark- ant iquark pair qq: 

|7>=|99>+- ( 2 - 47 ) 

The pair of qq forms a color dipole which scatters off the hadronic target. The dipole 
formulation of the rescattering problem is often called the s-channel approach to 
unitarisation which is alternative to the presented above t-channel approach. Here 
emission of a gluon can be viewed as a splitting of a dipole into two child dipoles 
in the s-channel producing a cascade of dipoles in s-channel. In the previously 
discussed EGLLA we introduced new degrees of freedom in t-channel and we required 



18 




the unitarity to be fulfilled in the t-channel and subchannels. In s-channel one starts 
with s-matrix S(itoi,b,Y) where x i=x — xi is the size of the dipole and x and 
xi are the transverse coordinates of the quark and antiquark and b = (xo + xi)/2 
is the impact parameter. We consider a frame where the dipole is right moving and 
the target is left moving and the all relative rapidity of the dipole and the target 
is taken by the target. This assures that there is only a small probability for the 
emission of a gluon from the dipole before the interaction. We may ask now how 
the S^Xoijfe, Y) evolve when we increase the rapidity by small amount dY. Here 
we decide to change (increase) the rapidity of the dipole which consequently emits 
a gluon. Due to the Fierz decomposition, the emission of a gluon in the large N c 
limit can be viewed as the splitting of a parent dipole into two daughter dipoles. 
The probability for the production of quark-antiquark-dipole state from the initial 
quark- antiquark dipole is given by: 

dP = ^ d x 2 dY-^ (2.48) 

When this probability is multiplied by the S matrix it the gives change in the S 
matrix so that one can write evolution equation [39]: 

dS(xoi*,Y) = aM I rf 2 X2 ^ [5(xQ2;b 1 jy)5(xi2jb _l 2 0j y)_ 5 (xoi,b,y)] 
oY 2ti z J XQ 2 xf 2 2 2 

(2.49) 
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where the first term corresponds to the scattering of the two dipole state on the 
target under the factorization assumption: 

(xoa, x 12 , b, Y) = S(x 02 , b + ix 12 , y)S'(x 12 , b - ix 20 , y). (2.50) 

This assumption is called mean field approximation to the gluonic field in the target. 
The exact treatment leads to an infinite set of coupled integro-differential equations 
where the equation for scattering of the two dipole state on the target involves the 
scattering of three dipoles. The other term corresponds to virtual contribution. 
One can rewrite that equation in more useful for our purposes form. Invoking the 
relation between the S matrix and the forward scattering amplitude N, S = 1 — N 
we obtain: 

dN(x 01 ,b,Y) a s N c f l2 xg, r , T/ , 1 _ „ T/ , 1 

QY = i^J d X2 ^[ 2 [N( * 02 ' b + 2 Xl2 ' Y) + iV(Xl2 ' b + 2 X20 ' Y) 

-iV(x 01 , b, Y) - iV(x 02 , b + ±x 21 , y)iV(x 12 , b + ix 20 , Y)\ 

(2.51) 

That formulation in terms of the forward scattering amplitude makes an easier 
link to approach in the i-channel. The direct comparison is, however, impossible 
because of different spaces of formulation. The impact factor in the momentum space 
corresponds to the Fourier transform of square of the wave function with the phase 
factors. The forward scattering amplitude N can be related via a Fourier transform 
to the dipole density in the momentum space which then can be further transformed 
to obtain the unintegrated gluon density. We discuss that in detail in chapter 3. 
The linear part of that equation is the dipole formulation of the BFKL equation [11]. 
Those terms are leading when the amplitude is small (i.e. the number of dipoles is 
small) and the quadratic term can be ignored. In the case of scattering off the nuclei 
one can interpret that contribution as single scattering of a single dipole off a nucleon 
in a nuclei. The nonlinear term introduces the possibility of simultaneous scattering 
of two dipoles off two nuclei and becomes important when the amplitude grows. That 
term contains the same TPV that we introduced before, but the most general proof 
(with dependence on momentum transfer) is not trivial due to complicated structure 
of the TPV. In the proof, Bartels, Lipatov and Vacca wrote most general fan diagram 
equation in configuration space. Then they performed the Fourier transform of the 
momentum expression of the TPV. The next step required restrictions on impact 
factors that they belong to the Mobius representation. With that assumption their 
fan diagram equation became equivalent to the BK equation [21]. We are going to 
discuss similar case later under the restriction that there is no momentum transfer. 
It has been shown [41, 59] that BK equation unitarises the cross section for fixed 
impact parameter and also suppress the diffusion into the infrared region which was 
draw bag of the BFKL equation. 
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It, however, violates the unitarity for a large impact parameter [41, 42] which is due 
to lack of a confining mechanism included in the BK equation. In this thesis, in part 
devoted to phenomenological applications we study properties of the BK equation in 
situations where the impact parameter dependence is provided by initial conditions. 
This allows us to avoid problems coming from diffusion to large dipoles. Further 
type of limitations are linked to the mean field approximation. The BK equation 
misses the effects of pomeron loops which are important at low densities. We show 
a how single pomeron loop can be constructed within the QCD, but the problem of 
introducing it into the evolution equation we leave as a task for the future. 
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Chapter 3 

The Triple Pomeron Vertex 



The Triple Pomeron Vertex in the perturbative QCD has attracted significant at- 
tention in the recent years. This particular is due to nonlinear evolution equations, 
e.g. the introduced previously BK equation, where the nonlinearity is given by the 
TPV. More recently, also the question of pomeron loops has been addressed: here 
again the TPV plays a central role. Whereas in many studies and applications it is 
convenient to use the coordinate representation, it is important to understand the 
structure also in momentum space. 

In this chapter we will investigate aspects of the TPV, related to its dependence 
on momenta arguments starting from the momentum space representation of the 
2^4 reggeized gluon transition vertex. So far, most of the studies of the nonlinear 
evolution equations have been done in the context of deep inelastic scattering where 
a virtual photon scatters off a single nucleon or off a nucleus. In both cases the 
momentum scale of the photon is much larger than the typical scale of the hadron 
or nucleus, i.e. one is dealing with asymmetric in configurations. As a first step in 
investigations the TPV, therefore, we will focus on the limit where the transverse 
momenta in the vertex are strongly ordered: this includes both the collinear and 
the anticollinear limits of the TPV. Furthermore we will also review and discuss the 
connection with other nonlinear evolution equation that have been discussed in the 
literature (GLR, GLR-MQ) 

3.1 The 2^4 gluon transition vertex 

The 2^4 gluon transition vertex introduced in previous chapter consists of three 
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Figure 3.1: Example of diagrams that contribute to the 2 — > 4 gluon transition vertex 
(wavy vertical lines represent reggeized gluons): real emission (left), a disconnected 
contribution (right). 



pieces: 



V^ saia3a9a *(Ki,K2;k a ,k 2 ,k3,k4) 



v/2~7T(5 a 'i a 2 

m - 1 



5 aia2 5 a3a4 y(K 1 ,K 2 ,k 1 ,k 2 ,k3,k 4 ) 



+5 a ^5 a ^V( Kl , k 2 , k 1; k 3 , k 2 , k 4 ) + 5 a ^5 a ^ V{k u k 2 , k 1; k 4 , k 2 , k 3 

(3.1) 

where Ki + k 2 = ki + k 2 + k 3 + k 4 = q and the subscripts a' { , a>i refer to the color de- 
grees of freedom of the reggeized gluons. It is convenient to express the 'basic vertex 
function' V(k,±, k 2 , k 1; k 2 , k 3 , k 4 ) with the help of the function G(k x , k 2 , k 1; k 2 , k 3 ) 



V(k x , k 2 , k x , k 2 , k 3 , k 4 ) = -g 4 



G(ki, k 2 , ki, k 2 + k 3 , k 4 ) + G(k u k 2 , k 2 , k x + k 3 , k 4 ) 



+G(ki, k 2 , ki, k 2 + k 4 , k 3 ) + G(k 1: k 2 , k 2 , ki + k 4 , k 3 ) - G(k u k 2 , ki + k 2 , k 3 , k 4 ) 
-G(ki, k 2 , ki + k 2 , k 4 , k 3 ) - G{k x , k 2 , ki, k 2 , k 3 + k 4 ) - G(k x , k 2 , k 2 , k l5 k 3 + k 4 ) 

+G(ki, k 2 , ki + k 2 , -, k 3 + k 4 

(3.2) 

The function G(k x , k 2 , ki, k 2 , k 3 ) [22, 23] generalizes the G function introduced in 
[6] to the non-forward direction. This function can again be split up into two pieces: 

G(k x , k 2 , ki, k 2 , k 3 ) = Gi(ki, k 2 , k x , k 2 , k 3 ) + G 2 (k x , k 2 , k x , k 2 , k 3 ), (3.3) 

where the first part describes the emission of a real gluon: 

r f i i i \ (k 2 + k 3 ) 2 ^ (ki + k 2 ) 2 K 2 \4k\kI 
Gi(ki, rc 2 , ki, k 2 , k 3 j = — , xn + 



v«i - ki 



k 2 - k 3 ) 2 (ki - ki) 2 (K 2 - k 3 ) 2 

-(ki + k 2 + k 3 ) 2 (3.4) 
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and the second one describes the virtual contribution: 

2 2 

G 2 ( KuK2 MMM) = -^(k(k 2 )-cu(k 2 + k 3 )]^ 2 )(^i-k 1 ) 

+ [u(k 2 ) - w(ki + k 2 )]5 (2) ( Kl - ki - k 2 )) (3.5) 

As it has been already stated the vertex in (3.1) is completely symmetric under the 
permutation of the four gluons. Moreover, it has been shown to be invariant under 
Mobius transformations [10], and it vanishes when Kj or kj goes to zero. 

This vertex can be used, for example, to construct a pomeron loop. A simple 
example for a pomeron loop in case of an elastic scattering of two virtual photons, 
is shown in Fig. 3.2(left). With the coupling of the pomeron loop to the external 
photons given by the impact factors, ai '° 2 , the expression for it reads: 

A(s, t) = -i2sn^ [ dY 3 [ dY 2 [ dY 1 5{Y -Y 1 -Y 2 - Y 3 ) 
2- Jo Jo Jo 

Swf^f^ 2 ^ q " ^^""hYs; K, Ki, q) 

/(£/ ^0^--(K 1> q-« 1 ;k 1> -k 1 -r > k il> -k 4 + r + q ) 

/ ^^ g(2)aia2blb2{Y2] kl ' k ' lJ *)<3 {2)aM ^ ks, k' 3 , r + q) 

/ S^ Vblb263b4;W(k - ^ ~ r ' k - " k 3 + r + * q ~ ^ 

g(2)mm( Yl . K ' v q)^(«', q - K '). (3.6) 

Here s is the squared center of mass energy, Y = ln(s/s ) is the total rapidity, Yi, 
Y 2 , Y 3 are the rapidity intervals as depicted in 3.2 which sum up to the total rapidity 
Y. Q^^M (2.20) is the (non-amputated) BFKL Green function which satisfies the 
BFKL equation. The combinatorial factor ^ in front of the integral is due to the 
symmetry under the interchange of the two pomerons which has its origin in the 
symmetry property of the 2 — > 4 vertex with respect to the interchange of outgoing 
gluons. The procedure that leads to (3.6) is outlined in the appendix. 

As a second example, we write the expression for the propagation of the BKP 
[8, 9] state between the the 2 — > 4 gluon vertices Fig. 3.2(right). Such a state 
contains the pairwise interaction of all four reggeized t-channel gluons. This is in 
contrast with the pomeron loop where gluons are forced to form pomerons therefore 
not all possible interactions are taken into account. Its Green's function satisfies the 
following evolution equation: 

(a; - w(k!) - u(k 2 ) - W (ka) - c(k 4 ))^ 4) {a!} ' K} ({k,}, {k^}) = 

0(4)0 KbK} ({kt} , {k ; } ) _ £ -±- 2 K^ 2 {b} ® G$ {hM<} ({h}{K}) (3.7) 

(ij) k ^ 
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Y, 



Y-2 



Figure 3.2: Contributions to the elastic scattering of virtual photons which contain 
the 2 — > 4 gluon vertex: a pomeron loop (left), a four gluon BKP state (right). Dark 
blobs represent Green's functions of reggeized gluons. 

where we have used the shorthand notation {kj} = (ki, k2, k3, etc. The sum 
extends over all pairs (ij) of gluons, and the kernel If^-iT includes the color tensor 
faibicfajbjc and the convolution symbol <g> stands for J The inhomogeneous term 
has the form: 

^ (2) (E k ^E k ^ (4) ° K}K} (^}'«}) = (2nrf[ 6a ^ {2) ^- K) . (3.8) 

1 K i 

The expression for the propagation of the BKP state between the 2 — > 4 vertices 
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reads: 



4 f Y f Y f Y 
A(s,t) = -t2sn- dY 3 dY 2 dY 1 5(Y — Yi — Y 2 — Y 3 ) 
4' Jo Jo Jo 



d 2 K <P Kl ^ _ K)g[ 2 Wl a' 2 ;a»a»(Y 3 ; K, Kl , q) 



/ (2tt)3 (2nY 

J n (§) ^we* - <o /n (H) (^)^(E^ - q) 

V" ; '* ai " 2a3 "Vi, q - Ki; kx, k 2 , kg, k 4 )^( 4 )^^H^} ( y 2 . {k,}{k;}) 
/ l^l^***^' k ' 2 ' k - k - ^ q - ^ 

£(2)6^^. ^ ^ qtyWtf, q - K >) (3.9) 

The statistical factor appears due to symmetry of the expression under interchange 
of the four gluons. It is easy to see that, when going from the BKP state to the 
state of two noninteracting BFKL Pomerons, the statistic factor ^ is replaced by 
I,. Again, a more detailed discussion of this diagram is presented in the appendix. 

In the following we focus on the pomeron loop and investigate, for zero momen- 
tum transfer q = 0, in the kinematic limit where, say, the momentum scale of the 
upper photon is much larger than the lower one. This implies that, at the upper 
TPV, the momentum from above, k±, is larger than the momenta from below, k 1; 
k3, or r ('collinear limit'). Conversely, for the lower TPV we have the opposite sit- 
uation: the momenta k^, k' 3 , and r are larger than k[ ('anticollinear limit'). Let us 
become a bit more formal. We expand the amplitude of Fig. 3.2 in powers of Q\jQ\ 
('twist expansion'). The object of our interest is the self-energy of the Pomeron 
Green's function, £(/Ci, k[). In eq.(3.6), £(ki, k[) is defined to be line 3-5, i.e. the 
convolution of the two TPV's with the two BFKL Green functions between them. 
It has the dimension k 2 , and it is convenient to define the dimensionless object 
£(§f) = S(Kl ' K p with the Mellin transform: 

poo 

E( 7 ) = / dtftik 2 )^ 2 )-*- 1 . (3.10) 
Jo 

The inverse Mellin transform reads: 

t{k2)= XS (A;2) ^ (7) ' (3 - n) 



where k = —h and the contour crosses the real axis between and 1 see Fig. 
3.3. Our analysis will then reduce to the study of the singularities of the function 
£(k 2 ). The twist expansion corresponds to the analysis of the poles located to the 
left of the contour in the 7 plane: the pole at 7 = is the leading twist pole, the 
pole at 7 = — 1 belongs to twist 4 and so on. For the upper TPV in Fig. 3.2, the 
analysis of the twist expansion requires the 'collinear limit', for the lower TVP the 
'anticollinear' one. 
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Figure 3.3: Singularities in the 7 plane. 



3.2 The collinear limit 



In this section we are going to study the collinear limit of the TPV. The ordering 
of the transverse momenta is the following: |k|^>|ki|, [k^j , | k3 1 , |k 4 | (and expansion 
parameters are |ki|/|k|, |k 2 |/|k|, |k 3 |/|k|, |k 4 |/|k|). In our investigations we will be 
interested in attaching color singlet objects to the vertex so we project (3.1) onto the 
color singlets using p a ^ a ^ =S aia2 5 asa4 /(N% - 1). In the multicolor limit iV c -> 00 
we obtain: 



pM v K;a }(Ki; k;j) = § a[,a' 2 V 2„ 



N* - I 

+V(1,4,2,3) 



(iV c 2 -l)y(l,2,3,4) + \/(l,3,2,4) 

^y(l,2,3,4) = vg'i(l ) 2 ) 3,4) 

(3.12) 



where ~5 a >* = V2n5<^ , V{1, 2, 3, 4) = V(k u k 2 ; k 1; k 2 , k 3 , k 4 ). 



3.2.1 The real part 

Let us begin the analysis by expanding the real part of the G function (3.3) in 
the collinear limit. In the analysis we are going to limit ourselves to the for- 
ward case. In the forward configuration we are going to use the simplified notation 
G(k, —k; ki, — ki — k 3 , k 3 ) = G(ki, k 3 ) and similarly for the vertex and basic vertex 
function V(k, — k; k 1; k 2 , k 3 , k 4 ) = V(ki, k 2 , k 3 , k 4 ). According to these definitions 
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the G\ (3.4) function reads: 

k 2 k 2 k 2 k 2 (ki + k 3 ) 2 k 4 

Gl(kl ' ks) = (k^k^ + (kfk^F " (k- kl ) 2 (k + k 3 ) 2 (3 - 13) 

where we used K\ = k, k 2 = — k. In the collinear limit the momenta of outgoing 
gluons in (3.13) satisfy conditions |ki|«|k|, |k 3 |«|k|. Performing the expansion 
up to second order terms we obtain from (3.13): 

G 1 (k 1 ,k 8 )~2k 2 



k r k 3 k r kk 2 2k 1 k 3 k 3 k 2ki-k3ki-k k 3 kk 2 
k 2 k 4 + k 4 k 4 + k 4 

2(k r k 3 ) 2 -k 2 k 2 



4 

(3.14) 

The first interesting term, which can be identified as the twist two contribution, is: 

G^k^-^k 2 ^ (3.15) 



where, since we are interested in convoluting the TPV with forward and angular av- 
eraged BFKL kernel we averaged over the azimuthal angle of k. Such contribution 
has the potential to give us the logarithmic integral over the transverse momentum. 
The presence of the momentum transfer would cause loss of a logarithmic contri- 
bution. We kept the same notation for the averaged function. After inserting this 
expression into (3.2) the whole expression vanishes. Passing to the twist four 



G 1 (k 1 ,k 3 )~2k 2 
we obtain for the TPV: 



2(k 1 -k 3 ) 2 -k 2 k 



Jl,2 



k 4 



(3.16) 



v r W (v v v \<\ jta ^ q ^ 4 u2 (- k i ■ k 2 k 3 ■ k 4 + kx ■ k 4 k 2 ■ k 3 + ki ■ k 3 k 2 ■ k 4 ) 

^LOAT c l kl ' k2 > k 3> k 4j —° 8 y k ^4 

(3.17) 

The expression above is the master formula for twist four contribution. In the 
above we introduced a shorthand notation. The superscript r stands for the real 
emission and similarly we will also use v for the disconnected (virtual) contributions 
in subsequent sections. In a spirit of our strategy we are going to consider the 
forward configuration of the TPV in the collinear limit i.e. when momenta are: 

ki = -k 2 , k 3 = -k 4 (3.18) 

Putting ki = 1, k 2 = —1, k 3 = m, k 4 = — m we obtain: 



28 



Figure 3.4: TPV with two interacting gluons attached to it(left). TPV with ladders 
of gluons (right). 




Figure 3.5: Kinematics used in 3.20. 

Having the expression for the vertex we are going to attach the BFKL ladders to 
its external legs. To achieve that, we multiply the vertex by propagators and act on 
it with the BFKL kernels (one for each leg). As a result we get the vertex and two 
interacting gluons attached to it. This procedure can be iterated and as a result we 
will obtain a vertex with attached BFKL ladders Fig. 3.4(right). We wish to obtain 
this result within leading logarithm accuracy and our goal is to get maximal power 
of the logarithm from the convolutions of the vertex with propagators and kernels. 
To do that we should act on the twist four contribution to the vertex with the twist 
four operator, which in our case is the product of two twist two contributions to the 
BFKL kernel (twist two for each leg). Let us compute the collinear approximation 
to the BFKL kernels which, convoluted with the TPV, will give the maximal power 
of the logarithm. The expression for the emission part of the BFKL kernel Fig. 
3.5 has already been presented, but we will recall it here (we will not write below 
K bfkl to simplify the notation): 

Assuming zero momentum transfer and q^»k^ we get: 

K coU ^-2kl (3.21) 

where we focused on the emission part since the virtual part does not contribute in 
the strong ordering approximation. The subscript coll stands for collinear. Using 
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this formula in the case of kinematics represented in Fig. 3.4(left) we get the fol- 
lowing expression for convolution of the vertex with the kernels (one kernel for each 
leg of the vertex): 



2<f 



d 2 \ 



d 2 m 2k 2 2k 2 , (2(l-m) 2 -l 2 m 2 ) 



2 hi (2vr)3 7 k , (2vr)3 l 4 m 4 



(3.22) 

where k is the lower scale which we do not specify at present. We notice that 
integrals over 1 and m give logarithms, but what is more striking is that angular 
integral over angle between 1 and m reduces the expression above and in particular 
the Triple Pomeron Vertex to zero. This result can also be obtained without ex- 
panding. To see that, it is enough to angular average over the azimuthal angles of 
momenta in the vertex (see appendix): 



V(k, -k;l,-l,m,-m)=4|- 



2k 2 #(l 2 - k 2 )6(m 2 - k 2 ) 



l 2 



+ In — 5(V - k 2 )6(m 2 - l 2 ) + In ===- 5(m 2 - k 2 )6(\ 2 - m 2 ) 



in 



r 



(3.23) 



The presence of 9 functions forbids collinear configuration. The physical meaning of 
that result is clear. If the two pomerons entering the vertex from below have smaller 
momenta than the pomeron from above, they cannot resolve it and cannot merge 
because they do not feel 'color' and the vertex vanishes 1 . 



3.2.2 The virtual part 

So far we have investigated contributions that come from the real part of the vertex. 
What is remaining is the disconnected part. Virtual pieces will be investigated 
convoluted with an impact factor. To deal with infrared finite quantities we choose 
to work with the impact factor of the photon. The function G 2 (ki,k 3 ) (3.5) in the 
forward direction reads: 

G 2 (k 1 , kg) = -2k 4 In |kl1 . 5(kl - k 2 ) - 2k 4 In -M—5(kl - k 2 ) (3.24) 

|Kl + k 3 | 1^1 + K 3 | 

and the photon impact factor (transverse which is the leading one) is: 

«W(k,Q) = S a '^a s a em ye 2 C drdp ^ - (1 ~ ~ (1 ~ ^ (3.25) 

In this expression p denotes the longitudinal component of k (in Sudakov decom- 
position) while second integration variable r is the Feynmann parameter. However, 



1 This has first been noticed in [34], see our discussion in chapter devoted to nonlinear evolution 
equations. 
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for our investigations it is enough to limit ourselves to twist expansion. To perform 
twist expansion of the impact factor one has to perform the Mellin transform of it 
with respect to k 2 /Q 2 . With the Mellin transform (3.10) one has: 



and we obtain: 



<t>a' ia > 2 {l) 



^ a , ia , 2 (k,Q) = j —i^—j aX2 ( 7 ) 



jq^ c r (-7 + 3) r(- 7 + 1) r( 7 ) r( 7 + 2) 
(-7 + 1) -7 + 1 7 r( 7 + §) 



(3.26) 



(3.27) 



Expanding that function and closing the contour of inverse Mellin transform to the 
left we obtain the following collinear expansion of the impact factor: 



<p aW2 (l,Q 1 ) = S< a '^,Qi)^^C 
while in the anticollinear case one has 



14 



T-3 ln (gf 



1 2 



1 2 2 / 1 2 ^ 2 



^(P, Qo) = <^0(P, Qo) * - jr (jg ) + 



14 4 /p 2 

~~9~~3 ln \Ql 



(3.29) 

where C = eja s a em . The momentum scales are ordered as follows \Q±\ ^> |1| 3> 
|p| ^ \Qo\- Qi stands for the virtuality of the photon which is larger then the 
momentum l 2 : while Qq stands for a virtuality which is smaller then p 2 . From now 
on we will use the notation: 0(k, Q) = 0(k). 

The virtual terms in the leading part of the vertex can be combined to give: 



W} (k) <g> V v ^ a N ](l, -1, m, -m) 



1 

(2tt) 2 2 



In 



l 2 



1 + m) 2 

2 



+ In 



l 2 



+ In 



m 



(1 - m) 2 

2 



(i,Qi) 



+ In 



m 



0(m,Qi) 
(3.30) 



2 / Y \2\ 



Assuming that |1| > |m| and inserting the twist four contribution to 0(1) (Cg(^p 
we can expand (3.30) in terms of |m|/|l|: 



M(k)®VlS(l,-l,m 1 -m) ~ -2C-^6^^ 



o m 2 l 4 (1-m) 2 l 4 
l 2 Q\ + l 4 Q\ 



. in- \ in 2 m 4 , / 1 • m \ L m 4 



(3.31) 
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After convolution with BFKL kernels we obtain: 

0K}(k) ® Vg^l, -1, m, -m) ® (K acolll K acoll2 ) ~ -2C^^^£J| In M 

1 (3.32) 

We see that this expression does not give us the expected power (third) of logarithm. 
Thus we arrive at the conclusion that we do not get a contribution from the TPV 
in the collinear limit in the leading logarithmic approximation. 



3.3 The anticollinear limit 



3.3.1 Real parts 

Let us now investigate the anticollinear limit. In the anticollinear configuration we 
allow for momentum to be transfered between the legs of the vertex. The momentum 
transfer here, as we will see, does not lead to loss of a logarithm. In this configuration 
w«|wi|, |w 2 |, |w 3 |, |w 4 | see Fig. 3.6. To study the emission part of the TPV it is 
convenient to rewrite the G\ function in the form: 



Gi(wi, w 3 ) = 2w 2 



1 



W-Wj 

'2 



W-W3 

'2 



^ _ w-wi -\2 



w-w 3 \2 



w-wi 

■w'f ■ 



W-W3 ' 



(3.33) 



the expansion parameters are |w|/|wi| and |w|/|w 3 |. Performing the expansion we 
obtain: 



Gi(w 1 , w 3 ) ~ 2w 2 



1 - 



w • W 3 W • Wi 



Wn 



w 



Wi-W 3 2 

2 2 W 
WJW3 



(3.34) 



Using (3.33) and averaging it over the azimuthal angle of w we obtain for the TPV: 



VlonAV, -P 



q, -q + r) ~ W^ 9 —1w 2 



(3.35) 



This term also appears when forward configuration is considered. Upon angular 
averaging over angles (without any expansion) this term survives, showing that the 
anticollinear configuration is preferred by the TPV (see appendix). However, in 
order to get the required logarithmic contribution after convolution with the BFKL 
kernels we need to consider higher order terms in 3.34. The resulting contribution 
is the following: 



V[ { ol(w 1 ,w 2 ,w 3 ,w 4 )^^^|-2w 4 



Wi'W 3 W 2 -W 3 Wi-W4 W2-W4 



w 2 w| 



w 2 w 2 



w 2 w| 



W2W4 



wi • (w 3 + w 4 ) w 2 -(w 3 + w 4 ) w 3 -(wi + w 2 ) w 4 -(wi+w 2 ) 



W 2 (W 3 + W 4 ) 2 W^(\V 3 + W 4 ) 2 W|(W! + W2) 2 W|(W!+W2) 2 

(wi + w 2 )-(w 3 + W 4 ) 
(Wi + w 2 ) 2 (w 3 + W 4 )^ 



(3.36) 
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Figure 3.6: Configuration of the momenta that is considered in the anticollinear 
limit. 



To proceed further we need the anticollinear limit of the BFKL kernel. Using (3.20), 
setting kx = Wi, k 2 = — Wj — r, q x = p, q 2 = q, — r and requiring that |q x | 3> 
|wi|, |r| we obtain: 

K acoU = -2w r w 2 , (3.37) 

where w 2 = — wi — r. As already mentioned before we are interested in a particular 
configuration, where there is a momentum flow between the legs of the vertex, but 
where there is no momentum flow out of the vertex. Such a configuration is realized 
by setting: W!=p, w 2 =— p— r,w 3 =q, w 4 =— q+r. In order to obtain the logarithmic 
contribution after the acting of the BFKL kernels on the TPV, we are going to 
consider the following cases: 



The configuration where |r|> |p| > |q|. Using (3.36) we obtain for the vertex: 

2 p 2 q 2 



V£(P, -P - r, q, -q + r) ~ -S^Aw^ (3.38) 



while the propagators read: 



1 1 



p 2 (p + r) 2 p 2 r 2 
1 1 



q 2 (q + r) 2 q 2 r 2 



(3.39) 
(3.40) 



using (3.37) we obtain for the kernels: 

K x = 2p(-p -r) ~ -2p r = K acolll (3.41) 

K 2 = 2q- (q + r) ~ 2qr = K acM2 (3.42) 
Convoluting those expressions we obtain: 

v^'}^^ u \ r ld2r f rf2 P f d ^ 42p-r2q-r2p.q 

V LONc ® (K acolll K acoll2 ) ^ 2 jJ w2 J^l 2 ^l 2 ^™ ^q^iv " 
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5 (2vr)9 2 3!\ ln ^J ' (3 ' 43) 

where |w | is specified only by the condition that it should be smaller than the 
momentum scales |1| and |m| which were considered in the collinear limit. Due 
to the additional integration over the momentum transfer this expression has 
different dimension than considered in the collinear limit. Convolution with 
the impact factor 0{ fe /}(w) = — b h '^\C In Ql/w 2 yields: 

5 L (27rP 2 4! \CQl 
(3.44) 

where in order to get the logarithmic contribution we took the lowest order 
term in the expansion of 0. In the configuration |q|>|p| the same result is 
obtained. 

Repeating the similar analysis in the case when |q|>|p|>|r|, we obtain: 

V> ® VS ® (^i^) = -^'^^pyf (l-f) 4 (3-45) 



And when |q| > |r| > |p|, we obtain: 



*W ® {KacMlKacMl) = ^Kjgg^ (l„ (3.46) 

The configuration |p| > |r| > |q| gives the same contribution. 



3.3.2 Virtual parts 

Let us now analyse the contribution coming from the virtual parts of the vertex 
in the anticollinear limit. Using (3.2) and (3.24) we list below a set of functions 
which contribute in the case when |p|, |q|> |r|: The other configurations will not 
contribute. 

<V} ® <? 2 (P, r) = -<^^p ^ |J^(P) " ln ^^ r ) ( 3 - 47 ) 



V) ® G 2 (-p - r > r ) = ln n^r 0((p + r)2) " sm ln p 0(r) 

(3.48) 
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G 2 (q, -r) = -^ 6: 



(27T) 



In 



q-r 



(q) - ^ fe 



' V 2. 



(27T) 



In 



q-r 



(3.49) 



/ } <g> G 2 (-q + r,-r) = -5 b i 



6', 6' 



(2^ 



In 



L0((q-r) 2 )-^ 



(2vr) ; 



In 



q 



(r) 
"(3.50) 



Analyzing the case when |q| >|p| >|r|, we obtain contributions from (3.47), (3.48). As 
an example, let us look for contribution arising from (3.47). Expanding 1/2 In ( - p ^ r - )2 </ ) ( r ) 
in terms of r 2 /p 2 we obtain: 



- In — ( 

2 (p + r) 2 



1 j - L 5 4r 2 



r 2 
p 2 



r P r 2 2 f r-p X 
p 2 p 2 \ p 2 



(3.51) 



Proceeding similarly with the other G 2 function we obtain 



{&'} 



V 



,r{6'} 
LON c 



,C2 i q 2 



2 5 (2tt) 2 r 2 



In 



(3.52) 



Convoluting it with kernels and setting the lowest scale equal to Qq we find: 



(3.53) 



which is exactly the same expression as (3.45), but with the opposite sign. The 
remaining G function gives contribution in the case where |p|>|q|>|r| which cancels 
againist with analogous coming from the real part. 

Let us shortly summarize our results before we continue the finite N c analysis. 
The contribution coming from the anticollinear part of the vertex is given by (3.44) 
when |r|>|p| >|q| and |r|>|q|>|p|. In the case when |p| > |q| > |r| and |q| > |p| > |r| 
we do not get contribution, because the appropriate terms cancel each other. If 
|q|>|r|>|p| and |p| >|r|>|q| we finally obtain (3.46). Our goal was to find the terms 
of the twist expansion of the TPV which, after convolution with the BFKL kernels 
would generate the maximal possible power of the logarithm. In the collinear case we 
expected to find two logarithms and after convolution with an impact factor, a third 
one. Disconnected pieces did not give contributions either. In the anticollinear case, 
however, we find the logarithmic contributions in the two of the specified regions. 
In the forthcoming section we are going to investigate contributions to the vertex 
that are suppressed in the large iV c limit. 
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3.4 Finite N c 



3.4.1 The collinear limit 

Analysis of previous sections repeated for subleading parts in N c gives: 
v r{a>} fn2 . 9 A 1 2 8(k r k 2 k 3 -k 4 - k r k 3 k 2 -k 4 + k^ka-ka) , , 

Substituting ki = 1, k 2 = —1, k 3 = m, k 4 = — m we obtain: 

, , ff 4 u2 81 2 m 2 
V s „ 6i ^ c (l, m, -1, -m) ~ ^ yk -j-j- (3.55) 

And analogously for the second subleading part. Analysis of the second and third 
non leading term give the same results so we limit ourselves to the second one. The 
convolution with the two BFKL kernels gives: 

(u u \^fW &± ,,29 A dH f* d 2 m 2kg 2kg (81 2 m 2 ) 

(2tt) 2 g 4 kgl / kV 
iV c 2 -l(27r)6 2 k 2 2! V ko / 
Convolution with the impact factor gives: 

r n ^ a i a i 2 C2tt1 3 rt 4 k 4 1 / (l 2 \ 3 

® ViU(^u^ co « 2 ) - ^ 5 I^F IT I 1 " W (3 - 57) 

As we have already observed in order to obtain the logarithmic contribution from 
disconnected piece, we are forced to expand G function in terms of ratio |l|/]m| or 
reverse. It turns out that the virtual pieces of subleading in N c part of the vertex 
give contribution. The contribution comes from the following functions: 



ha'} ® G 2 (l, -1 - m) = -5«-±-\n JlU(l) - 5 a '^-±- 2 In + m) 

(27T) 2 |m| (2tt) 2 |m| 

(3.58) 

K} ®G 2 (-l,l+m) = -^ ia2 ^ ln ^#)-^' iai ^ ln ^^^( 1 + m ) (3-59) 

® G 2 (m, -1 - m) = -<^_I_ i n ^0(m) - <^^p In + m ) 

(3.60) 

K} ® G 2 (m, 1 + m) = -^' ia ' 2 ^p In - <^^p In ^p0O + m ) 

(3.61) 
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Expanding the functions above in terms of |m|/|l|, collecting all terms and integrat- 
ing over angles yields for convolution of Vj^^ c (l, m, —1, — m) with the impact factor 
and kernels, the following result: 

0K> ® V S U (W^) * -—(JC--^^ ^ (ki^ j (3.62) 

The power of the logarithm is the same as in the real part, when we compare (3.63) 
with (3.57). The total contribution reads: 

**> 8 V -= 8 - -ife^Tgfll ( to §) 3 (3 ' 63) 

3.4.2 The anticollinear limit 

Here the results differs only in the kinematical region where |q|>|p|>|r| and |p|> 
|q|>|r| i.e. we do not get a piece from the vertex that convoluted with kernels give 
logarithms. The other regions give the similar results as the leading N c . To be 
explicit: 

• |r| > |p| > |q| and |r| > |q| > |p| 

0M ® ® (K acM1 K acM2 ) * -^^^^f (in 4 (3.64) 

which is exactly the same as (3.43), but suppressed by N% — 1. 

• |p| >|r|>|q| and |q|>|r|> |p| 

*m ® vS » = -j^jIj^ y f ("> j) * (»■«) 

which is exactly the same as (3.46), but suppressed by JV| — 1. 

Repeating the similar analysis for the virtual part we observe that none of the 
specified regions gives the leading logarithmic contribution and the real part gives 
the total contribution. The considered limits do not cover all possible cases. We 
could for example consider situations when the momentum transfer is the largest 
scale in the problem. 



37 



Chapter 4 



Nonlinear evolution equations 



4.1 The Balitsky-Kovchegov equation for the un- 
integrated gluon density 

In this section we are going to derive the BK equation. Let us consider the state of k 
reggeized gluons in the Heinsenberg picture which is labeled by color and momentum 
degrees of freedom: 

\k) = |ki..k fc ;ai,..a fc ) (4.1) 

The normalization is: 

k 



1 

! 

cr(fc) 1=1 



(k\k>) = S Wk - J2 II (( 2 ^( k * - K)^) (4-2) 



where the sum is over the permutations of outgoing gluons. Using this states as a 
base we assume now that the proton state can be written as: 

oo 

fc'=l 

Evolution of this state is given by: 

e- yH \ P ) = \ P (y)), (4.4) 
where the Hamiltonian is given by: 

H = H 2 ^ 2 + H 2 ^ 4 + H 4 ^ 2 (4.5) 
Its matrix element are defined via: 

(k\H 2 ^ 2 \k') = 5 kk , [t^^t E E fa^cfc^K^ k v K K) fl 5 ^ - k '^ aib ' 1 

'' a(k'-2)i>j=l l^ij 

k' k< 



i I 

(4.6) 
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This part corresponds to the BKP interaction. The last two terms account for gluons 
which propagate freely without interaction. 



k k' 

(k\H^ 2 \k') = 5 kk ,„ 2 J2 J>-W(k, 1 k*;k; ) ^k| ) k:) 

s>t i>j>t>r=l ^ 

5(k^ + k; + k; + k;-k s -k t ) f[ 5(k p -k' p )5 a ^ 

this part of the Hamiltonian corresponds to the changing number of gluons interac- 
tion. It allows for four gluon to fuse into two gluons. 

Finally the part which corresponds to the transition of two gluons into four gluons 
is given by: 

k' k 

(k\H 2 ^\k') = 5 kk , +2 £ V aia i aia ^<(k h kj, k h k r ; k' s , k' t ) 

s>t=l i>j>t>r=l . . 

k> ^ 

Siki + kj + kt + K-ki-kj) 11 5(k p -k' p )5 a ^ 

Let us define the wave function of k gluons in the proton at rapidity y in the following 

way: 

jfi} = {k \e-y H \p) = (2 7 r) 3 ^ } (y,k 1 ,k 2 ...k n ) (4.9) 
Upon differentiating it with respect to y we obtain: 

= -(k\He- yH \p) = -^{k\H\k'){k\e- yH \p) 

y k' 

k' 

The unity operator is given by: 

£*WI = £II / ^^a^al (4.11) 



k< i=l 

where the summation on the left hand-side means also integration over the continu- 
ous degrees of freedom. The expression (4.10) is a set of the infinite, many coupled 
equations. It cannot be closed because, for instance equation for two gluon wave 
function involves contribution coming from four gluon wave function. 

r\<paici2 

= -((2\H 2 ^ 2 \2) ® J=)\^ - «2|tf 4 _ 2 |4> ® F,) a ^ (4.12) 

The term proportional to (4.8) vanishes since it requires zero gluons in the initial 
state. The minus sign appearing in front of the linear term is due to our definition 
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of the BFKL kernel. In forward case one gets the usual relative sign difference 
between linear and nonlinear parts. One can, however, write closed equation using 
factorisation ansatz: 

^ 4 ( 2 /,k 1 ,k 2 ,k 3 ,k 4 ) aia2a3a4 = c^ 2 ( 2 /,k 1 ,k 2 ) aia2 ^ 2 ( 2 /,k3,k 4 ) a3a4 (4.13) 

This assumption is strictly speaking justified for a large nuclei where it is natural to 
consider coupling of one pomeron to one nucleon. At present we, however, assume 
that we can use this approximation for the proton therefore we obtain: 

= -((2\H 2 ^ 2 \2) <g> J= 2 ) a ^ - c((2|#4^ 2 |4> ® {T 2 T 2 )) aia \ (4.14) 

The factorisation (4.13) gives some freedom of choosing the parameter c [32]. We 
will choose the constant c equal to l/y/2 to get the usual version of BK equation. 
To obtain the BK equation for the unintegrated gluon density let us define the 
unintegrated gluon density via: 

rr 2 {vMM) = (2 7 r) 3 ^ 2 (y,k 1 ,k 2 )^- = (k\e-y H \p) (4.15) 

In the next step we act on (4.14) with 5 aia2 (such color tensor is present for instance 
in photon impact factor). Using (4.6), (4.7), (4.8), (4.14), (4.15) we obtain: 

dT 2 {x, k,q) f d 2 \ F 2 (x,\,q) 



/d 1 J~2\X 
77r-^K BFKL (\,ci- l;k,q- k)-^— 
{2ny r(q-l) 
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(4.16) 



1 2 (-1 + § + r) 2 m 2 (-m + f - r) 2 



Factors Cp arising from contractions of color tensors were absorbed in definition of 
function JF 2 . The factorial indicates the symmetry of the TPV. We put ki = k and 
for the second momentum argument k 2 =— ki + q we wrote just q to indicate that 
there is a dependence on the second momentum variable. The momenta incoming 
into the vertex from below is labeled k' : = 1, k' 2 = — 1 — q/2, k' 3 = m, k^ = — m — 
q/2 — r. The variable r stands for the loop momentum. If we restrict ourselves to 
the zero momentum transfer at the top of the fan diagrams (which corresponds to 
neglect the impact parameter dependence) and, inside the fan diagrams, carry out 
the integrations over momentum transfer under the factorization assumption: 

F 2 (x,k,r)=F(x,\t.)F(r,R), (4.17) 

where in the present case 

F(r,R) = e -^— (4.18) 



40 




Figure 4.1: Diagrammatic representation of the integral version of the equation 
(121). 



F(r, R) is the form factor of the proton. The parameter R has meaning of a radius 
of the proton. After assumptions done above we obtain somewhat simpler equation: 



dJ r {x 1 k) 
d In 1 / x 
f d 2 l d 2 m 
2!"2tt# J (2tt) 3 (2tt) 3 



-7T- 



V(k, 



Kbfkl(i, -i; k, —k) — ^ — 
T{x, 1) T{x, m) 



(2tt) 3 
-k; 1, 



(4.19) 
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The assumed factorization limits the momentum transfer across the ladders to small 
values. Note that this procedure has been developed for deep inelastic scattering 
on a large nucleus and that nonzero momentum transfer allows to keep track of the 
impact parameter dependence in configuration space. In the next step we perform 
integrations over azimuthal angle of 1, m and k. We use f(x, k 2 ) for angular averaged 
function J-(x, k). After this procedure the linear part of the equation (4.19) is just 
the BFKL equation for the unintegrated gluon density while in the nonlinear part 
the basic vertex function takes the form (see appendix for details of this calculation): 
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(4.20) 



The fan diagram equation reads: 
df(x,k 2 ) iV c a Sl2 



d In 1/x 



71 



dK, 2 



K 



at 



f(x,K 2 )-f(x,k 2 ) , f(x,k 



+ 



(4K 4 + k 4 )J 
50 dm 2 



2R 2 



2k 2 / -rf x,l 2 / —f z,m 
A 2 1 A 2 m 4 

2 
2 



+f(x,k 2 ) ^ r ln^J/(x,l 2 ) + /(x,k 2 ) 



dm 2 



m 1 



In 



m l ft 2 



(4.21) 



41 



When applying this equation to the scattering of a virtual photon on a nucleus 
one may ask what is the most dominant contribution. In the DGLAP approach (see 
chapter 2) one has strong ordering in momentum, i.e virtualities of gluons closer to 
photon quark vertex are larger. Here (on a level of the evolution equation) we can 
also consider such configuration where the transverse momenta are strongly ordered: 
the upper momenta should be smaller than the lower ones. However, making use 
of our results for the collinear limit of the TVP and of the structure of the angular 
averaged vertex, where the 9 functions does not allow the collinear configuration, 
we arrive at the conclusion that, there is no such contribution. In contrast, the 
momenta are ordered in the opposite direction, i.e. at the nonlinear vertex of the 
equation, the momentum above have to be smaller than those below. In more phys- 
ical terms, the recombination of two smaller gluons ends up in a larger gluon. Note 
that this only holds, after the large- N c limit has been taken and the angular aver- 
aging has been performed. 



4.2 Comparison with other equations 

One can obtain (4.21) from the BK equation under the assumption that the target is 
large in comparison to the size of the dipole which scatters on it. The crucial point 
is the space of formulation and identification of the correct degrees of freedom. The 
BK equation was formulated in configuration space so the first step to compare it 
with (4.21) is to perform its Fourier transform. The next point is to notice that the 
equation obtained via Fourier transform has to be further transformed in order to 
take form of (4.21). This is achieved by invoking relations linking dipole density, 
dipole cross section and the unintegrated gluon density. To derive this one has 
to apply relations between a scattering amplitude of a color dipole 7V(x 01 , b, x), 
the dipole cross section <t(x i,x) (where x i is a size of the initial dipole) and the 
unintegrated gluon distribution f(x, k 2 ) [43, 44, 27]. Let us start with the Balitsky- 
Kovchegov equation for the dipole amplitude. 

N f x 2 1 1 

\ d 2 x 2 ^\- [iV(x 02 , b+-,x) + iV(x 12 , b + -x 20 , x) 

T J X 02 X 12 ^ ^ 

-iV(x 01 , b, x) - iV(x 02 , b + ^x 21 , x)iV(xi2, b + ix 20 , x)\ 

(4.22) 



<9iV(x 01 , b, x) 
d In 1 jx 
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The assumption that target is large in comparison to the size of dipoles allows us to 
simplify BK equation and rewrite it in slightly different form: 



jV(x 01 ,b,x) = a s N c r ^ 
dlnl/x 7r 2 J p 
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d 2 x 2 ^^iV(x 02 , 6, x) 7V(x 12 , 6, x) 
x 02 x 12 



(4.23) 



what allows to write the dipole amplitude in a factorised form: 

7V(x,b,x) = N(x,x)S(b) 



(4.24) 



where 5(b) is the profile function of the proton related to form factor (4.18) via 
two dimensional Fourier transform. The factorization ansatz will be used later. 
Performing Fourier transform of (4.23) [14] we obtain [14, 33]: 
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(4.25) 

(4.26) 



01 



d 2 k , 



and inverse transform is given by: 

iV(x 01 , b, Y) = x 2 x J ^-: e ik - X01 $(x, k, b) 
Multiplying both sides of (4.25) by k 2 and using Mellin representation: 

PCX) 

7 , b) = / rfk 2 (k 2 )-^ 1 k 2 $(x, k 2 , b) 
Jo 

we obtain: 

<9$(x, 7, b 



(4.27) 



(4.28) 
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where: 
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[2^(1) -#y)- 7(1-^)] 



(4.30) 



is the BFKL eigenvalue. To find the equation for the unintegrated gluon density we 
apply relation between the dependent on the impact parameter gluon density and 
the dipole scattering amplitude: 

N(xoi,b,s) = -JL J ^[1- J (|k||x 01 |)]aJ(x,k 2 ,b) (4.31) 
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where 



J d 2 bf(x, k 2 , b) = J d 2 bf(x, k 2 )S(b) = f(x, k 2 ) (4.32) 
using (4.26) we obtain: 

$(aab) = \ J ^^J^ J^ [1 - Mk^)]a s f(xX,h) . (4.33) 



1 f d 2 x i ?1 . b 8tt 2 f dk 

Integration over angles yields: 

^2 r°° a\A r°° 



o 2 roo 7i 2 /.oo j 

$(a;,l 2 ,b) = — / " / ^J (|l|)[l-J (|k||x 01 |)K/(a;,k 2 ,b), (4.34) 
and the integral over x 01 gives 

7t 2 r°° dv 2 (v 2 \ 

1 2 > b) = - jf a — In (^j «./(*, k 2 , b) . (4.35) 



Now we need to invert the operator 

/-oo ,71,2 



° = — /„ ir'"(p> (k) ' (4 36) 

where g(k 2 ) is a test function. Multiplying both sides of (4.37) by l 2 and performing 
the Mellin transform with respect to l 2 we obtain [27]: 

$(x,7,b) = J d\ 2 \ 2 <$>{x, l 2 , b)^ 2 )-^ 1 = 

j^iE*^)^**^ 1 ^'^ 1 <4 37> 



(1-7) 2 ' 



and equivalently 



N. 



f{x, 7 , b) = — £.(1 - 7 ) 2 $(a:, 7, b) . (4.38) 



The inverse Mellin transform gives the momentum space version of this expression: 
, 2 UA ^ c /• d 7 2 2 iV c 2 d , 2l 2. 



/(x,l 2 ,b) = -±j / -^(1^(1 -7)^(7, = -^(l-l 2 -2) 2 l^(x,l,b) . (4.39) 
a s ix z J 2m a s ii z dl 

Inserting (4.37) in (4.29), and inverting the Mellin transform we obtain: 

W^) = ^<3/(x,k 2 ,b) (440) 
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This equation can be further transformed, applying the factorization ansatz f(x, k 2 , b) 
f(x, k 2 )S'(b) and taking for example a Gaussian profile function: 



5(b) 



e r 1 



(4.41) 



After using this and integrating over the impact parameter, we obtain: 
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At this stage we would like to comment on the relation between (4.42) and 
the other known nonlinear evolution equations appearing in considerations of high 
energy limit of QCD. The first nonlinear evolution equation which was a milestone 
in physics of saturation is the GLR-MQ [34], [36] (Gribov, Levin, Ryskin, Mueller, 
Qiu) equation, which is given by (equation number (2.41) in [34] and (30) in [36]): 
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(4.43) 



<91n(l/a;)<91nk 2 ' 

in the equation above the constant C is not the same in both of cited papers, however, 
for our discussion its value is not crucial. This equation can be rewritten for the 
unintegrated gluon density with the collinear approximated BFKL kernel: 



N c a s 



dfjx, k 2 ) 
d In 1/x ix 



dk 



12 



■f(x,k' 2 )-C 



k' 2 

K o 



'2 



k 2 R 2 



dk 



f(x,k' 2 ) 



k 



(4.44) 



The nonlinearity is claimed to be given by the TPV at the collinear limit. From our 
analysis one, however, comes to the conclusion that at collinear limit (in angular 
averaged case) TPV vanishes and only anticollinear pole contributes. This can also 
be read off from comparing nonlinear terms in (4.42) and (4.44). The structure of 
integrals is totally different and one cannot transform one into another. The equation 
that has similar properties to (4.42) is the GLR equation (equation number (2.108) 
in [34]) derived from summing up at low x (single logs are summed) fan diagrams. 
It is written directly for the unintegrated gluon density (the gluon density in the 
GLR notation includes propagator). That equation is an attempt to generalize the 
BFKL equation for physics of dense systems. GLR obtained TPV (they did not 
consider disconnected pieces, the pre factors are different and they work in large N c 
from the beginning) and approximated it by a constant. This equation reads: 
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where g is the TPV vertex in the local approximation. The relation to integrated 
gluon distribution is given by: 



xg(x,Q 2 ) = 




dk 2 &(x,k 2 ) 



(4.46) 



Gribov, Levin and Ryskin before applying local approximation to their nonlinear 
part of their equation wrote expression for the TPV convoluted with gluon densities 
from below: 



The structure of this expression (where to have direct comparison one should replace 



formulation) is similar to the structure of the real emission contributions in the 
nonlinear part of (4.42). The lack of the disconnected contribution does not lead to 
singularities since (4.47) is infrared safe. From this structure, GLR arrive at similar 
conclusion as we that in order for gluons to fuse they have to resolve individual 
gluons in the ladder from above. 

It will be interesting to compare precisely all of the listed above equations and 
to check validity of approximations which were done. At this moment one can 
guess that although very different structure of nonlinearity in (4.42) and (4.44) the 
behavior of solutions will be similar in the saturated phase because the nonlinear 
terms have similar scaling. More precise studies are, however, necessary. 



a s (m 2 )a s (\ 2 )&(x, m 2 )&(x, l 2 )fl(l 2 -k 2 )fl(m 2 -k 2 ) 

(4.47) 
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)/k 2 , we did not do that since we want to refer to the original 
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Chapter 5 



Phenomenological applications of 
the fan diagram equation 

5.1 Model of subleading corrections to BK 

For phenomenological applications of BK it is desirable to have a formalism which 
embodies the resummation of the subleading corrections in In 1/x. Attempts in this 
direction already exist, see for example [43, 44, 57, 59, 63, 64, 65]. Equation (4.42) 
contains the BFKL kernel at leading logarithmic (LOlnl/x) accuracy. This is a 
coarse approximation as far as a description of the HERA data is concerned. It 
is well known [40] that the NLOlnl/x corrections to the BFKL equation are quite 
large. To make the equation more realistic, one can [43, 44, 27] implement in the 
linear term of (4.42) a unified BFKL-DGLAP framework developed in [2]. This 
framework unifies the BFKL and the DGLAP equations in alternative to a CCFM 
way [62] . This unification is however, on a different level of rigor then the CCFM 
(direct summation of diagrams). Kwiecihski, Martin and Stasto observed that the 
BFKL equation may be refined if one includes relevant parts of the DGLAP splitting 
functions and other contributions of subleading in In 1/x order. This contribution 
leads to coupled integral equations which at low x limit become equivalent to the 
BFKL and at collinear limit become equivalent to DGLAP evolution for the gluon 
density. Let us be more precise. In the KMS scheme [2], the BFKL kernel becomes 
modified by the consistency constraint Fig. 5.1 [47, 48] 

k /2 < k 2 /z , (5.1) 

imposed onto the real-emission part of the kernel in Eq. (4.42) 

f^dk' 2 r /(f,k' 2 ,b)fl(^-k' 2 ) - /(f,k 2 ,b) /(f,k 2 ,b) | 

./,, k' 2 1 k' 2 k 2 + |4k' 4 + k 4 |U- 

The consistency constraint on the real gluon emission(5.1) resumes a large part of 
the subleading corrections in In 1/x [45, 46]. The physical meaning of it is that the 
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virtuality of the gluon is dominated by its transverse momentum. 



x, k 




kx — k' T 



Figure 5.1: Kinematic variables appearing in the kinematic constraint 

Additionally, the non-singular (in x) part of the leading order (LO) DGLAP 
splitting function is included into the evolution: 



—K®f^> —K®f + dzP gg (z)f(-,k' 2 ,b) 

Jx Z Jx Z J K J X z 



where 
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2N r 



(5.3) 



(5.4) 



Additionally, we assume that in our evolution equation a s runs with scale k 2 which 
is another source of important NLOlnl/x corrections. The improved nonlinear equa- 
tion for the unintegrated gluon density is as follows [44]: 
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(5.5) 



above we assumed the target to have cylindrical profile with radius R = 4GeV 1 

8(R 2 - b 2 ) 



S(b) 



nR 2 



(5.6) 
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and we integrated over the impact parameter. In the KMS approach the inhomoge- 
neous term was defined in terms of the integrated gluon distribution 

/(°)(,,k 2 ) = ^Lj\ zPM lg (^X) (5-7) 

taken at scale kg = lGeV 2 . This scale was also used as a cutoff in the linear 
version of the evolution equation (5.5). In the linear case this provided a very good 
description of F 2 data with a minimal number of physically motivated parameters, 
see [2]. The initial integrated density at scale k^ was parameterized as 

xg(xX) = N(l-xY , (5.8) 

where N = 1.57 and p = 2.5. 

Let us finally note that we include most of subleading corrections only in the 
linear part of the BK equation. The nonlinear term gets contribution only from 
running of the coupling constant. Effects of energy conservation in nonlinear term 
are important as has been shown in analysis of BK in the configuration space (but 
subleading) [66] 

One could try to implement the kinematical constraint effects in the TPV in a 
similar manner as in the linear part i.e: 



V <g> f(x, k 



2kV: r - " ; 



This contribution is not considered here and we postpone it to further studies. 



5.2 Numerical analysis 

5.2.1 The unintegrated and integrated gluon densities 

In this section we introduce the method of solving Eq. (5.5) and we present the 
numerical results for the unintegrated gluon distribution function f(x, k 2 ) and the 
integrated gluon density xg(x,Q 2 ). The method of solving (5.5), developed in [44], 
relies on reducing it to an effective evolution equation in Inl/x with the boundary 
condition at some moderately small value of x (i.e. x = x ~ 0.01). 
To be specific, we make the following approximations: 

I. The consistency constraint 6(k 2 /z — k' 2 ) in the BFKL kernel is replaced by 
the following effective (z independent) term 

9(k 2 /z - k' 2 ) - 6(k 2 - k /2 ) + (Jjpj 6(k' 2 - k 2 ) . (5.9) 
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This is motivated by the structure of the consistency constraint in the moment 
space, i.e. 



J^ z -e(kyz-k' 2 )^e(k 2 -^) + (^j e(k ,2 -k 2 ), (5.10) 



2. The splitting function is approximated in the following way 

dz 



jf ™[zP gg (z) - 2N c ]f g,k' 2 ) - P gg (u, = 0)f(x,k' 2 ) , (5.11) 

where P gg {ui) is a moment function 

r 1 dz 

W= / -z"[zP gg (z)-2N c ], (5.12) 
Jo z 

and _ 

^> = 0) = --. (5.13) 

This approximation corresponds to retaining only the leading term in the 
expansion of P gg {oj) around uj — 0, see [49]. 

Using these approximations in (5.5) we obtain [27]: 
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(5.14) 

First, the equation (5.14) was solved with the non-linear term neglected starting from 
the initial conditions at x — 10~ 2 given by (5.8). The parameter u> c s was adjusted 
in such a way that the solution of the linear part of (5.14) matched the solution 
of the original equation in the BFKL/DGLAP framework [2]. This procedure gives 
^cfr = 0.2 and the solution of the linear part of (5.14) reproduces the original results 
of [2] within 3% accuracy in the region 10" 2 > x > 10" 8 and 2GeV < k 2 < 10 6 GeV. 
This matching procedure has also the advantage that the quark contribution present 
in the original BFKL/DGLAP framework is effectively included by fitting the value 
of uj e &. The full non-linear equation (5.14) was then solved using the same initial 
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conditions and setting R = AGeV" 1 . In Fig. 5.2, we plot the unintegrated gluon 
distribution function as a function of x for different values of k 2 . This figure com- 
pares the results of two calculations, based on the linear and nonlinear equations. 
The differences are not significant, however there is some suppression due to the 
nonlinearity at the smallest values of x < 1CT 5 . The subleading corrections strongly 
decrease the value of the intercept with respect to the LOlnl/x value and the non- 
linear term becomes important only at very low values of x. As is evident from 
Fig. 5.3, the subleading corrections cause a large suppression in the normalization, 
also at moderate values of x. 
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1 e+03 1 




Figure 5.2: The unintegrated gluon distribution f(x, k 2 ) obtained from Eq.(5.14) as a function of 
x for different values k 2 = 10 2 GeV 2 and k 2 = 10 3 GeV 2 (up) and for k 2 = 1.1 GeV 2 , k 2 = 5 GeV 2 , 
A; 2 = 50 GeV 2 (down). The solid lines correspond to the solution of the nonlinear equation (5.14) 
whereas the dashed lines correspond to the linear BFKL/DGLAP term in (5.14) 
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Figure 5.3: The same as Fig. 5. 2 (for clearer presentation we skip the plot with k 2 = l.lGeV 2 ) 
but now the modified BK equation (5.14) (solid lines) is compared with the original BK equation 
(5.14) without subleading corrections (dashed lines). 
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Figure 5.4: The unintegrated gluon distribution f(x,k 2 ) as a function of k 2 for two values 
of x = 10~ 5 and 10~ 4 . Up: solid lines correspond to the solution of the nonlinear equation 
(5.14) whereas dashed lines correspond to linear BFKL/DGLAP term in (5.14). Down: solid lines 
correspond to the solution of the nonlinear equation (5.14) whereas dashed lines correspond to the 
solution of the original BK equation without the NLLl/x modifications in the linear part (5.14). 
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Figure 5.5: The integrated gluon distribution xg(x,k 2 ) as a function of x for values of Q 2 = 
10 2 GeV 2 and Q 2 = 10 3 (up) and for Q 2 = 5 GeV 2 and Q 2 = 50 (down) obtained from integrating 
f(x, /c 2 )Eq.(5.14). Dashed lines correspond to solution of linear BFKL/DGLAP evolution equation. 
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This is due to the fact that the non-singular in x part of the P gg splitting function 
was included into the evolution. This term is negative and is important at large 
and moderate values of x. The same conclusions can be reached by investigating 
the plots in Fig. 5.4 where the unintegrated density is shown as a function of the 
transverse momentum k 2 for fixed values of x. The nonlinear effects seem to have a 
moderate impact in that region. On the other hand the subleading corrections are 
substantial. For example, at x — 1CT 5 and k 2 = 10 GeV 2 the reduction in magnitude 
of the unintegrated gluon density is about 25. Fig. 5.5 shows the integrated gluon 
density xg(x,Q 2 ) = dk. 2 f(x, k 2 )/k 2 . The change from the power behavior at 
small x is clearly visible in the nonlinear case. Additionally the differences between 
the distributions in the linear and nonlinear case seem to be more pronounced for 
the quantity xg(x,Q 2 ). This is due to the fact that in order to obtain the gluon 
density xg(x, Q 2 ) one needs to integrate over scales up to Q 2 including small values 
of k 2 , where the suppression due to the nonlinear term is bigger. 



5.2.2 The saturation scale Q s {x) 

In order to quantify the strength of the nonlinear term, one introduces the saturation 
scale Q s (x). It divides the space in (x, k 2 ) into regions of the dilute and dense 
partonic system. In the case when k 2 <Q 2 s (x), the solution of the nonlinear BK 
equation exhibits the geometric scaling. This means that it is dependent only on 
one variable N(r,x) = N(rQ s (x)) or in momentum space f{k 2 /Q s (x)). To calculate 
the saturation scale in our model with resumed NLlnl/x effects we follow [65]. There, 
the saturation scale was calculated from the formula 



This formula has been obtained by interpreting the quantity <&(x, k 2 ) or equivalently 
f(x, k 2 ) as field describing particle with coordinates (ml/a;, lnk 2 /ko). With such 
a field one can accociate group velocity and phase velocity. Studying the solution 
of linearized BK one realizes that in the particular case, when phase velocity and 
wave velocity become equal, the equation admits so called wave front solution. This 
happens for some anomalous dimension 7 C (7 is an argument of Mellin transform of 
the kernel). It was realized that the slightly modified traveling wave solution of the 
linear equation can be used as ansatz for nonlinear equation. The crucial point is 
that the exponent entering this ansatz is j c . It has interpretation of critical exponent 
relevant for saturation. Using this exponent it is possible to calculate effective slope 
of pomeron and saturation scale. For LOlnl/x BFKL one obtains 7 C = 0.23 The 
effective Pomeron intercept uj s is a solution to the equation 

Ws(l) = a a x(l,Us) , (5-16) 
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Figure 5.6: Saturation line obtained from (5.16, 5.17). 

where x{li u ) * s the kernel eigenvalue of the resumed model. In our case the eigen- 
value has the following form 

X(7,w) = 2tt(l)-¥( 7 )-tt(l-7 + W ) + ^P ro (a;). (5.17) 

The solution for the saturation scale obtained from solving (5.15,5.16) using eigen- 
value (5.17) is shown in Fig. 5.6 and gives A = = 0.30,0.45,0.54 for three 
values of a s = 0.1,0.2,0.3, respectively [27]. These results are similar to results 
obtained in [65]. We compare our results with the saturation scale from the Golec- 
Biernat and Wiisthoff model. Normalization of the saturation scale is set to match 
GBW saturation scale at xq = 0.41 x 10~ 4 . The saturation scale Q s (x) can be also 
obtained directly from the numerical solution to the nonlinear equation by locating, 
for example, the maximum of the unintegrated gluon density as a function of rapid- 
ity. At large x for fixed gluon momentum BFKL/DGLAP effects lead to a strong 
growth of the gluon density. At certain value of x the nonlinear term becomes equal 
to linear and cancels out with it. That effect leads to occurrence of a maximum 
which leads to a natural definition of a saturation scale. We are going to discuss 
this issue further in context of the gluon density with explicit dependence on im- 
pact parameter. For the purpose of phenomenology we attempt here to estimate 
the effect of the nonlinearity in a different, probably more quantitative way [27]. 
We study the relative difference between the solutions to the linear and nonlinear 
equations 

\f^(x,Q a (x,P) 2 )-f^(x,Q a (x,!3) 2 )\ 
f^(x,Q s (x,f3y) 

where (3 is a constant of order 0.1 — 0.5. Since this definition of the saturation 
scale is different from the one used in the literature and is likely to posses different x 
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dependence, we denote it as Q s . In Fig. 5.7(up) we show a set Q s which are solutions 
to Eq. (5.18) for different choices of (3 together with the saturation scale calculated 
from the original saturation model by Golec-Biernat and Wusthoff [50]. Solid lines 
given by Eq.(5.18) show where the nonlinear solution for the unintegrated gluon 
starts to deviate from the linear one by 10%, 20%,. . . , 50%. It is interesting that 
contours Q s (x) defined in (5.18) have much stronger x dependence than saturation 
scale Q s (x) defined by Eq.(5.15) and the one from GBW model. In particular 
Q s (x,/3) > Q s (x) for given x (at very small values of x). This might be a hint that 
saturation corrections can become important much earlier (i.e. for lower energies) 
than it would be expected from the usual definition of the saturation scale Q s (x). In 
Fig. 5.7(down) we also show contours in the case of the integrated gluon distribution 
function, that is the solution to (5.18) with f(x,k 2 ) replaced by xg(x,Q 2 ). As 
already seen from the previous plot, Fig. 5.5, the differences in the integrated gluon 
are more pronounced. For example in the case of Q 2 = 25GeV 2 and x ~ 10~ 5 — 10~ 6 
we expect about 15% to 30% difference in the normalization. Again, by looking 
solely at the position of the critical line, one would expect the nonlinear effects to 
be completely negligible in this region since at x — 10~ 6 the corresponding Q 2 {x) ~ 
2.8GeV 2 (taking Q 2 (x) = Q 2 (x/xo)~ x with normalization Q 2 = lGeV 2 at x ~ 
4 x 10 -5 and A ~ 0.28, [50]). This rough analysis shows that one cannot think 
of saturation scale as a definite and sharp border between very dilute and dense 
system. The transition between these two regimes appears to be rather smooth 
and the nonlinear term of the equation seems to have quite a large impact on the 
normalization even in the 'linear' regime defined as Q 2 3> Q 2 (x). In practice, the 
estimate of the saturation effects is even more complicated since the unintegrated 
gluon density has to be convoluted with some impact factor, and the integration 
over the range of scales must be performed. 
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Figure 5.7: Solid lines show contours of constant relative difference between solutions to linear 
and nonlinear equations, Eq. (5.18). Lines from bottom to top correspond to 10%, 20%, 30%, 
40%, 50% difference. Left: Contours in the case of the unintegrated gluon distribution f(x,k 2 ); 
right: contours in the case of the integrated gluon distribution xg(x, Q 2 ). Dashed line in both case 
corresponds to the saturation scale from Golec-Biernat and Wiisthoff model [50] . 
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5.2.3 Dipole cross section a(r, x) 

It is interesting to see what is the behavior of the dipole cross section a(r, x) 
[27] obtained from the unintegrated gluon density via: 

a(r,x) = —j ^(1- J (|k||r|))a/(:r,k 2 ), (5.19) 

this equation follows from (4.37) via integration over impact parameter. In this 
calculation we assume that a s is running with the scale k 2 and we use r size of 
the initial dipole, i.e. r = x i. Calculation of the dipole cross section requires the 
knowledge of the unintegrated gluon density for all scales < k 2 < oo. Since in 
our formulation the unintegrated gluon density is known for k 2 > k 2 we need to 
parametrise f(x, k 2 ) for lowest values of k 2 <kg. We use the matching condition 

/•k 2 J k 2 

xg(xX) = J ^/(^k 2 ), (5.20) 

and following [51] we assume that f(x, k 2 ) ~ k 4 /kQ for low k 2 . This gives (compare 
Eq.(5.8) ) 

f(x,k 2 )=AN(l-xY(k/k ) 4 . (5.21) 

In Fig. 5.8 we present the dipole cross section as a function of the dipole size r for 
three values of x = 10~ 3 10~ 4 and 10~ 5 . For comparison we also present the dipole 
cross section obtained from GBW parametrisation. To be self-consistent, we cut 
the plot at r = 2 GeV~ x because we assumed in the derivation of formula (4.42) 
that the dipoles are small in comparison to the target size (we assume proton radius 
to be 4GeV _1 ). This cut allows us to obtain a model independent result since we 
observe that different parametrisation of f(x, k 2 ) for k 2 < kg give essentially the 
same contribution for r < 2GeV _1 . We observe that our extraction of the dipole 
cross section gives similar result to the GBW parametrisation. The small difference 
in the normalization is probably due to the different values of x g which probe the 
gluon distribution (or alternatively the dipole cross section). In the GBW model 
the dipole cross section is taken at the value x g = x which is the standard Bjorken 
x = Q 2 /2p ■ q. On the other hand, in the formalism presented in Ref. [2] one takes 
into account the exact kinematics (energy conservation) in the photon impact factor. 
It is a part of the subleading effect in the impact factor and it increases the value of 
x g ~ 5x. Therefore, in our formalism the normalization of the unintegrated gluon is 
increased so that the convolution with the impact factor and the resulting structure 
function remains the same. 

5.3 Description of the F 2 structure function 

In this section we are going to use the BK equation with subleading corrections to 
describe the F 2 structure function [67]. To do so, first of all we are going to use 
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Figure 5.8: The dipole cross section obtained from modified BK (solid line) compared 
to GBW dipole model (dashed line). 



more realistic profile function of the nucleon. Instead of the cylindrical profile of the 
nucleon we are going to use Gaussian one. Our motivation is to control the impact 
parameter dependence in a more detail in order to have more general formalism 
applicable for prediction for more exclusive processes as for example diffractive Higgs 
production. The description of F2 is performed to fix parameters of the input. The 
Gaussian profile is more realistic than the cylindrical since it introduces dependence 
of screening corrections on the distance from the center of the target and does not 
lead to over-saturation of the gluon density at very large rapidity. With that setup 
the BK equation reads: 

df(x,k 2 ,b) 



N c a s (k 2 ) k2 

7T 



dk 



12 



k 2 
K o 



dhv(l/x) 
f(x,k 2 ,b) (6 (k 2 - k' 2 ) + v k , : 



k 2 \ ^eff 



£(k' 2 -k 2 ))-/(*,k 2 ,&) 



f{x,k 2 ,b) 



[4k ,4 + k 4 ]i 



|k /2 -k 2 | 



a.Jk 



ZTT J 1.2 



dk- 

Tv2 



f{x,k' 2 ,b) 



-2vra s (k 2 ) 



2k 2 



r°° dk 2 



f(x,k' 2 ,b) +2f(x 



, 2 r, r dk 



12 /1J2' 



f(x,k 2 ,b) 

(5.22) 



the inhomogeneous term stands for the input gluon distribution and is given by: 

f i0 \x,k 2 ,b) = ^S(b) f dzP ga (z) X -g (^,k 2 ) . (5.23) 
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Figure 5.9: The unintegrated gluon density as a function of b for fixed k 2 and x 
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As usual the unintegrated gluon density is obtained via the integration over distance 
from the center of the target: 



f(x,k 2 ) = d 2 bf(x,k 2 ,b) 



the input profile is 



S(b) 



exp( 



-b 2 - 



7T R 2 



(5.24) 



(5.25) 



We take R = 2.8GeV~ 1 which follows from the measurement of diffractive J/i/j 
photo-production off proton. In fig (5.9) we plot the unintegrated gluon density 
as a function of the impact parameter for fixed values of k 2 and for x = 10~ 8 . We 
observe that at small values of x and central values of b the saturation effects are 
strongest and lead to the depletion of the gluon density. On the other hand, due to 
large distance phenomena as for example confinement which we model via Gaussian 
input, there are less gluons for peripheral b. The net result is clear. The impact 
parameter dependent gluon density has a maximum as a function of distance from 
the target. 

Fig. 5.10(down) visualizes the unintegrated gluon density as a function of x. At 
large x for fixed gluon momentum BFKL/DGLAP effects lead to a strong growth of 
the gluon density. At certain value of x the nonlinear term becomes equal to linear 
and cancels out with it. This effect leads to occurrence of a maximum which leads 
to a definition of a saturation scale [67] i.e. we define the saturation scale as Q 2 for 
which: ^j^r^ = . Using that definition of saturation line we plot critical line 
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Figure 5.10: The unintegrated gluon density as a function of x for b=0 and fixed k 2 (up), The 
unintegrated gluon density as a function of x for 6 = 3 and fixed k 2 (down). 

for three values of impact parameter (5.11). The extracted value of parameter A 
according to prescription Q 2 S = (f) X gives A = 0.293 for b 2 = OGeV- 2 , A = 0.298 
for b 2 = 3GeV~ 2 and A = 0.298 for b 2 = 6GeV~ 2 . Those values agree with the 
A parameter obtained in previous section where we studied saturation line coming 
from wave front analysis. It was obtained for a s — 0.1 and A = 0.29. That value 
is also very close to the one in GBW model which we already referred to. It is 
worth to note that the normalization that we obtained in case of b 2 = 0GeV~ 2 was 
not far from the GBW model i.e. Xq = 0.45 x 10~ 4 while that value in GBW was 



63 




100 F 



k 2 =1.2GeV 2 



calculated with Gaussian profile 
calculated with cylindrical profile 



10 



M 




°le-08 le-07 le-06 le-05 le-04 le-03 le-02 

x 

Figure 5.11: Critical line obtained for three values of the impact parameter (up). The uninte- 
grated gluon density integrated over impact parameter for different profile functions. Continuous 
line corresponds to Gaussian profile while scatter line corresponds to cylindrical target. 

0.41 X 10" 4 . 

A similar maximum is not seen for the gluon density integrated over the impact 
parameter Fig. 5. 11 (down), it flattens, but does not fall as very small x is ap- 
proached. This is due to the large contribution to the integral from the peripheral 
region where density of gluons has not saturated yet. In Fig. 5.10 we observe well 
known fact that the lower k is the earlier saturation effects manifest themselves. 
This can be understood in our approach from the structure of the integral in the 
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nonlinear part. The lower limit of integration is given by value of k 2 at which the 
gluon is probed and extends to the infinity. The lower the momentum is, the longer 
is the path of integration. Let us now proceed to calculate the F 2 structure function. 
The method of calculation is the following. We use the solution of (5.22) to obtain 
the unintegrated gluon density which depends on the distance from the center of 
the target. Then we integrate it over the distance from the center of the target. 
The next step is application of the KMS prescription to calculate F 2 (x, Q 2 ). In this 
approach, i.e. the sea quark contribution is calculated within the kr factorisation 
including NLOlnl/x corrections coming from exact kinematics. The valence quarks 
are taken from [70]. Within the KMS prescription: F% ea = Yl q $q : where: 

where the impact factors corresponding to the quark box contributions to gluon- 
boson fusion process are the same as those used in ref. [2] (see also [7]), i.e.: 

S b q ox (z, k, Q) = a em [ d(3d 2 K'a s 5(z - z ) (5.27) 
Jo 



{[3 2 + (1 - Pf 



K n — k 



D 



lq 



D 



2q 



+ [m +4Q^(l-/3) 



where k' — k — (1 — (3)k and 



D 



lq 



D 



2q, 



(5.28) 



D lq = k 1 + (3{1 - (3)Q 2 + m\ 



D 



lq 



{K-k) 2 + P(l-P)Q 2 + m 2 q 
k 2 



n' 2 + m 2 



-i 



*= l+ w^w + wA (5 - 29) 

The effect of kinematical NLOlnl/x corrections is included via argument of unin- 
tegrated gluon density in formula (5.26) which is defined as x g = x/z. This is an 
effect coming from exact kinematics of the gluon-boson fusion. At leading order 
we would have z = 1 in (5.29) and the integral over z would (neglecting quark 
masses) reduce (5.29) to leading LOlnl/x formula (2.30). So, at leading order we 
have (f> q (k,Q) = f z dz^S q ox (z,\s. 2 ,Q 2 ), where 0(k, Q) was introduced in the first 
chapter. 

The full NLOlnl/x corrections are still not known, but there is activity in that di- 
rection [68, 69]. In the calculation of the (effective) quark distributions we use the 
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Figure 5.12: Diagrammatic representation of the F 2 calculation. 



impact factors (5.28) corresponding to the massless quarks and the (charmed) quark 
mass effects are included in the threshold factor: 



a c (k 2 ) = 1 + 



k 2 + m 2 



(5.30) 



In the impact factors we use (5.28) with m u = = m s = and m c = lAGeV. 
The various components entering the box in the KMS prescription are treated in 
the following way 

• region when k 2 , k,' 2 < kg is assumed to be dominated by the soft pomeron 



exchange 
where 



q<j qP I qP i qP 

D u ' d ' s 



= S% = 2Sf = 0.037 ln(462/a;) (1 - x) 



(5.31) 



(5.32) 



where we modified the power like contribution from nonperturbative region in 
order to be consistent with saturation picture 
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• k 2 < kg < k' 2 the strongly ordered approximation is applied, what results in 
neglecting k 2 in comparison to Q 2 and n 2 

• The region k 2 > kg is treated perturbatively and perturbative expression for 
S g is used 




Figure 5.13: F 2 structure function. 

The charm component of £ is generated perturbatively, i.e. there is no soft pomeron 
contribution to it and it is given by: 

S q=c (x, Q\ )= - S box (z, k 2 = 0, Q 2 ; m c )^-/(-, k 2 ) 
Jx z Jo k z 



+ 1 — S box (z,k 2 = 0,Q 2 ;m c )-^/(-,k 2 ) (5.33) 

z Jkl K Z 



With that formalism we obtained good description of data which made it possible 
to fix parameter of the input gluon distribution: 

xg{x, ko) = 1.572(1 - x) 2 5 (5.34) 
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The effect of nonlinearity is, however, hardly visible and is just at the edge of HERA 
kinematical region Fig. 5.13. We also would like to stress that the effect of saturation 
is better visible and more pronounced in case of the impact parameter dependent 
gluon density. This observation suggests that nonlinear effect may easier easy to 
observe for quantities which are sensitive to that parameter or its Fourier conjugate: 
the momentum transfer. 

5.4 Rescattering corrections in the diffractive Higgs 
boson production. 

In this section we will compute the hard rescattering correction to diffractive central 
Higgs production at LHC energies [71]. The diffractive Higgs production in collision 
of protons is of interest since the background is much smaller then the expected 
signal. The other advantage is that the mass of produced Higgs boson can be 
measured in two independent ways; first with the high accuracy by measuring the 
missing mass to the forward out going protons and, second by the H — > bb decay. 




Figure 5.14: Kinematics of the Two Pomeron fusion contribution to the exclusive 
Higgs boson production. 
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Standard approach 

In the standard pQCD [72]- [78] approach, it is without including the hard rescat- 
tering corrections the amplitude for process under studies is given by: 

M (y) = i2^A J ^J° ff {xu k 2 ; »)f off (x 2 , k 2 ; p) (5.35) 

Here, the constant above is given by: 

A 2 = KVG F a 2 s (M 2 H )/9ir 2 (5.36) 

where Gp is the Fermi coupling constant, and K m 1.5 is the NLO K— factor. The 
quantity f°^(xi, k 2 ; fi) is so called off-diagonal gluon distribution and is given by: 

d 



dQ 2 



xg(x,Q 2 )^T g (Q,fx) 



(5.37) 
Q 2 =k 2 



where T g (k,fi) is the Sudakov form-factor and it accounts for probability that there 
will be no multiple emission from fusing gluons, 

_2 2 ^ r( A + 5/2) 
^ " /? T(A + 4) (5 ' 38) 



is the correction factor and allows us to obtain off-diagonal gluon distribution if we 
now diagonal, the \x scale is // ~ 0.62M# ■ The last important point is the so called 
opacity factor Q(s, b) which accounts for nonperturbative soft rescattering. Those 
effects are factorised from the hard process leading to the Higgs production under 
assumption that the hard process leading to the production of Higgs occurs on short 
time scale. Assuming this factorisation is fulfilled the effect of soft physics enters 
just by multiplication of the hard amplitude by so called gap survival factor which 
depends on the opacity S 2 = exp(Q(s,b)). The differential cross section following 
from the amplitude (5.35) (without including gap survival factor) is: 

da pp ^ pH p(y) = \M (y)\ 2 
dydt x dt 2 256tt 3 1 ' ' 

After inclusion of the proton elastic form factors one obtains the dependence of the 
cross section on momentum transfers t\ = (pi — ki) 2 and t 2 = (p 2 — k 2 ) 2 of the form 
of exp(i? 2 (ti + 1 2 )/2). Thus, after the integration over momentum transfers t\ and 
t 2 , the single differential cross section reads: 

da PP ^ P Hp(y) = \M (y)\ 2 

dy 647T 3 /? 4 ' 1 ' 

For our future developments it is convenient to have the right-hand side of expres- 
sion above in impact parameter representation using the Gaussian profile (5.25) we 
obtain: 

^ (0 %Wp = _1_ j d2b j d 2 hils(hi)s(h _ bl )M (y)| 2 , (5.41) 
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where b x and b 2 = b — hi are transverse positions of the production vertex measured 
from centers of the protons and b is the impact parameter of the collision. 

5.4.1 Hard rescattering corrections 

The motivation to include hard rescattering corrections to diffractive Higgs produc- 
tion comes from the large rapidity gap between the colliding protons. At rapidity 
19Y apart from single logarithms of Higgs boson mass small x logarithms should 
be resumed. To be specific let us consider process where two protons scatter at 
large rapidity. According to our previous discussion of saturation interacting par- 
tons may overlap and recombine what will have impact on the Higgs production. 
That effect may be sizable because of large rapidity available. We would like to 
stress here that we are going to consider only a subclass of the screening corrections 
which should be included in BFKL-based description of the Higgs production. In 
particular we are going to restrict ourselves to pomeron fan diagrams configurations 
in the exchanged system which leads to Balitsky-Kovchegov equation for the gluon 
density. The considered configuration is therefore highly asymmetric and should be 
considered as estimate of importance of rescattering corrections. To be specific the 
amplitude consists of following elements: 

• The two scale unintegrated gluon distribution, f( a \xi, ki; /i a ) describing glu- 
ons above the TPV. This gluon distribution is evaluated at relatively large X\. 
The scale /i a = ki 

• The two-gluon Green's function Ga0<-2, k 3 ; %2', l^a-, l^b) that describes the 
propagation of the two gluon singlet state between the Higgs production vertex 
and the TPV. The diffusion approximation for the Green's function is going 
to be considered GbfklQz-2, k 3 ; xi, X2', /if,)- The Green's functions will have 
typical fib ~ ran /2, which belongs to the Higgs production vertex. In order 
to account for the scale dependence at the Higgs vertex, the Sudakov form 
factor, y/T g {kk \ /if,) will be included into the Green's function: 



where D = U((3) 

• The Higgs boson production vertex which does not differ from the TPF 

• The off-diaginal unintegrated gluon distribution, f b (x3, k 3 , q — k3; gluon 
distribution that enters the Higgs boson production vertex from the lower 
side, with the total transverse momentum q. That gluon distribution is taken 
to be the same as in the case of the TPF i.e. it is an off-diagonal gluon 
distribution given with /if, = 



k 2 ||k 3 | y/T g (k 3 ,/if,) Mog^kj/kjn /xi 

2vr v /nDlog(x 1 /x 2 ) \_ADlog(xi/x2)\ \x 2 



) 



A 



Gbfkl0<-2, k 3 ; xi, x 2 ; /if,) 



(5.42) 
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Figure 5.15: The hard rescattering correction to the exclusive Higgs production. 



• The off-diagonal unintegrated gluon distribution, p c '{x±, k 4 , — q — k 4 ; /i c ), for 
the screening Pomeron that propagates across the large rapidity distance from 
the lower proton to the TPV with momentum transfer — q. The scale \x c is set 
by the gluon virtuality, so one has f <yC \xA 1 k 4 , — q — k 4 ; k 4 ). 

After these definitions the dominant part of the diagram can be written: 

dx4 f d 2 q f d 2 ~k 2 f d 2 k 3 f d 2 k 4 



2 J ^""^3 j 27rk 4 
— -\f {a \xi, ki) V(ki, -ki, k 2 , -k 2 , k 4 , -k 4 )G A (k 2 , k 3 ; ari, x 2 \ m H /2) 

f^(x 3 , k 3 , q - k 3 ; m H /2)f^(x 4} k 4 , -q - k 4 ) (5.43) 

where: 

A/f A/f k^ k^ k^ 2 

x 3 = -^-e~ y , x 2 = -^=-e~ y , xi = ^ and x a ~ -j-, x& ~ The scale k that enters 
the definitions Xj corresponds to the average values of the virtualities, we set that 
scale k^ = lGeV 2 . The scattering is assumed to be forward, however, due to the 
fact that below the TPV the loop is formed the momentum transfer (q) dependence 
enters / (a) . The dia gram accounting for replacement of target an projectail is 
accounted by replacing y — > —y. The correction amplitude is then: 

M corr (v) = M^riv) + M«.(-y) (5.44) 

The integral above is infrared safe due to emergence of the saturation scale Q 2 (x 4 ) in 
gluon f( c \ That scale at interesting us values of x 4 as HERA data indicates is well 
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above the cut-off scale kg. To implement the soft rescattering correction which will 
enter through so called opacity factor it is needed to have formula for the differential 
cross section in impact parameter representation. It reads: 

da(0+1){ dy PP ^ PHP = -^ J d% J d 2 h 1 \S(h 1 )S(b-b 1 )M (y) + M cmr (y,h,h 1 )\ 2 , 

where the dominant parts of both M (y) and 

d 2 q f d 2 k 2 f d 2 k 3 f d 2 k 



(5.45) 



ImMW r (y,bM) = -IGnVA H ^ [ 

Jx a X 4 J 



(2tt) 2 J 2^k\) 27r\4J 2nk\ 



I ^Tl[/ (a) (^i, b - bi) <g> V](k 2 , k 4 )G A (k 2 , k 3 ; x u x 2 ; m H /2) 
J 27rk 1 

/W(x 3 ,k3,b 1 ;m H /2)/( c )(a; 4 ,k 4 ,b 1 ) - {y ^ -y}(5.46) 

The impact parameter dependence of and is given by S(b) and may be 
factorised out, in analogy with the case of the Two Pomeron Fusion amplitude. The 
impact parameter dependence of is determined by the BK evolution, and it is 
different from 5(b) because of nonlinearity. Having defined M corr (y, b, bi), we can 
write down the cross section for the exclusive production, taking into account the 
soft rescattering 

da {0+1) ' n (v) 1 r r 

pp ^ HpVy) = — J d 2 hj d 2 b, |5(b 1 )5(b-b 1 )M (y)+M COIT (y,b,b 1 )| 2 exp(-Q(. 

V 71 (5.47) 

where Q(s, b) is the opacity of the pp scattering at squared energy s and impact 
parameter b. Here we adopt the opacity factor of the two channel eiconal model 
proposed in [79, 80]. 

5.4.2 Hard rescattering corrections - results and discussion. 

All numerical calculations of the amplitude and of the cross section for the central 
exclusive Higgs boson production have been performed choosing y = and assuming 
the LHC energy y/s = lATeV. The mass of the Higgs boson was assumed to be 
M H = 120GeV. In Fig. 5.16 we show ratio of corrected amplitude to TPF amplitude 
as function of transverse distance bi of the Higgs boson production vertex from the 
center of the target proton: 



R b (h) 



M corr {y,b,bi) 
M (y,b,b 1 ) 



(5.48) 

y=0,b=0 



The amplitudes are integrated over y 1 . We observe that the ratio (5.48) decreases 
with increasing distance, bi from the proton center. Thus, the corrections term 
tends to suppress the Higgs production inside the target proton. The correction 
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Figure 5.16: The ratio of the hard rescattering correction to the Two Pomeron Fusion 
at b—0, as a function of the transverse, distance, b\ of the Triple Pomeron Vertex 
from the center of the target proton, for the two values a s = 0.15 and a s = 0.25 at 
^7s = UTeV. 



g?o"tpf / dy 


da /dy (a s = 0.15) 


da /dy (0.17) 


da/dy (0.2) 
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is large for both choices of a s . In Table 1 we show values of differential da/dy 
cross section for TPF and TPF with rescattering corrections included. Results for 
different values of a s are presented. At a s = 0.2 the absolute value of the hard 
rescattering correction exceeds the Two Pomeron Fusion term and this explains 
why, at this value of a s , the corrected cross section takes its smallest value. For 
all considered values the cross section with rescattering corrections is smaller then 
TPF and is sensitive to the triple pomeron coupling. As we have already said we 
considered only subclass of diagrams that contribute to the described process. The 
more sophisticated considerations would require fan diagrams from both sides [81] 
and the inclusion of the pomeron loop or BKP states. 
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Chapter 6 

Summary and conclusions 



In this thesis we have studied aspects of the high energy limit of QCD. We investi- 
gated theoretical problems related to the physics of saturation and its implications 
for phenomenology. 

In the second chapter we constructed amplitudes for elastic scattering of hadronic 
objects with the exchange of a pomeron loop and also four gluon BKP state. To find 
preferred by the TPV regions of phase space we studied collinear and anticollinear 
limits of the Triple Pomeron Vertex. This analysis allowed us to gain some insight 
into the complicated structure of the vertex. In particular using this analysis one 
can see which kind of scattering proceses are mediated via triple interaction. Our 
analysis has shown that in the collinear limit the TPV does not give the expected 
power of logarithm when convoluted with BFKL ladders and impact factors. We 
find this result striking and the possible explanation is that the preferred configu- 
ration for the system is the configuration where there is no triple interaction when 
the scales are strongly ordered. The preferred configuration in that situation may 
be when gluons directly couple to the quark loop. We also studied the anticollinear 
limit of the TPV and we found that the leading contribution does not give the ex- 
pected power of logarithms and comes from the anticollinear pole. The situation 
changes at higher orders and also when we allow momentum to be transfered then 
there is a contribution with the expected power of logarithm. Similar analysis was 
carried out also for a finite number of colors. We show that in this case there is the 
expected power of logarithms both in the collinear and the anticollinear limit. 
The result involving infinite number af colors was also obtained by performing angu- 
lar angular averaging over the azimuthal angles of momenta variables in the vertex. 
We found that the real emission part consists of a term which corresponds to the 
anticollinear pole accompanied by theta functions indicating anticollinear ordering 
of scales. The virtual term (which is subleading) also includes a theta function 
indicating anticollinear ordering of scales . This result is in agreement with the 
performed twist analysis. In chapter four we applied the Triple Pomeron Vertex 
together with the BFKL kernel to propose the Hamiltonian for the evolution of the 



74 



proton which we assumed to be state of reggeized gluons. We derived coupled evo- 
lution equations for the system of these gluons. These equations are very similar to 
those studied by Balitsky. We however, propose them in momentum space where 
the interpretation in terms of Feynmann diagrams is much clearer. We show that 
assuming mean field approximation which basically means large N c limit we obtain 
nonlinear closed equation which is equivalent to the the BK equation for the unin- 
tegrated gluon density. Finally we discuss the relation of BK for the unintegrated 
gluon density to other known nonlinear equations. 

The fifth chapter is devoted to phenomenology. We develop a procedure to make 
the gluon density that follows from BK more realistic in order to apply it to phe- 
nomenology. The important point is to find an efficient way to implement subleading 
corrections in Inl/x corrections and to choose realistic proton profile function. We 
follow the method of KMS which is easy to implement and is consistent with exact 
NLOlnl/x calculations. Included corrections have a large impact on the behavior 
of the gluon density. They are especially important for large value of momentum. 
The effect of nonlinearity is not so important at this region. This is due to the fact 
that the path of integration in the nonlinear term is short. It however, becomes 
important for smaller values of the gluon momentum and small x where it leads 
to the depletion of the gluon density (for cylindrical target). Similar behavior is 
not observed for a gluon density with a Gaussian profile (integrated over impact 
parameter) since at largest distances the effects of saturation are weaker. For that 
reason we argue that the Gaussian profile is more physical since it does not lead 
to over-saturation of the gluon density. With the introduced setup we perform a 
fit to F 2 data to constrain gluonic input distribution. The form of the input gluon 
distribution can be still chosen to be more optimal. At present we choose it to be 
power like and all depletion is due to the nonlinearity of the equation. We also 
study the saturation line that follows from BK and we show that the saturation 
effects on level of the gluon density are almost negligible in the HERA kinematical 
region. We wish to check in the future analysis another functional form of it. Finally 
we develop the method to include hard rescattering corrections to Diffractive Higgs 
boson production at the LHC. We argue that at the rapidity available at LHC one 
can expect hard rescattering corrections to be important (this follows from studying 
the saturation line). Those corrections should be treated separately from soft ones 
which were already studied. The crucial points here are: 

• the presence of the TPV which allow pomerons to rescatter 

• constrained by HERA data the gluon density coming from the solution of 
(5.22). 

The obtained corrections turn out to be rather large and decrease the cross section 
for Higgs production. We however, keep in mind that those results are a rather 
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rough estimate. To be more realistic one should consider more symmetric situation 
where fan diagrams are allowed to propagate from two sources. One should also 
include in this processes contributions coming from closed pomeron loops (3.6) or 
even BKP states(3.9). The task to include these diagrams in their full complexity 
is difficult. We can however, try to investigate in a future analysis the importance 
of the contribution from the different parts of phase space and apply our results of 
collinear analysis to simplify the problem. 
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Appendices 



A Details of the pomeron loop calculation 

In this appendix we will give short outline of sequence of steps that led us to (3.6, 
3.9). The general structure is transparent and reflects the Regge factorization. To 
set up the convention for the impact factors let us begin with the quark-quark 
scattering due to the color singlet exchange. The averaged over helicities and colors 
amplitude in eiconal approximation equals to: 



4(1) 
qq^qq 



>,t = 0) = Usn J ^C^C (A-l) 



where as usual \ti is the momentum of t-channel gluon and s is the square of the 
total energy, (j) a q ia2 = S aia2 ^-. The next step is to calculate the elastic scattering of 
quarks at three loop level due to the color singlet exchange, Fig. 6.1 (left). Including 
s and w-channel contributions we arrive at: 

„m / x ,„ „ 1 , s f d 2 kx d 2 k 2 d 2 k 3 1 1 1 

NIK bfkl (\l u \l 2 )K bfkl (\l 2 , k 3 )4>° m 

where ki, k2, k3 are momenta of t-channel gluons. 

To investigate the formation of the pomeron loop it is enough to calculate two 
loop correction to the process above. The correction, that we are interested in, comes 
as two additional gluons in the t-channel exchanged between gluons in the s-channel, 
Fig. 6.1 (right). We are in particular interested in color singlet contribution which 
therefore does not contribute to reggezation. In calculation of that contribution we 
make use of dispersion technics. The result is: 
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Figure 6.1: Left: amplitude for quark quark scattering due to color singlet exchange 
in high energy limit QCD (A. 2). Right: central cut given by given by (A. 4). 



Disc M A^ qq (s,t= 0) = 2 S 7r 2 27rP^^ fc ' 1 ^P a 3^^^_±_i.l n 2 ± 

< ra\a2 

f d 2 k, d 2 k 2 d 2 k 3 d 2 r ^ijjg^^ 

J (2vr) 3 (2tt) 3 (2vr) 3 (2tt) 3 (2tt) 3 k^(k 2 -r) 2 (k' 2 + r) 2 ^ 

0a" ia " 2 (A.3) 

The combinatorial factor (^y) 2 in the formula above reflects symmetry with respect 

to interchange of gluons in left and right- hand- side parts of considered diagram. 

The factor ^ comes together with logarithm of the energy and indicates the process 

of the exponentiation. Other contributions as peripheral cut may be obtained by 

means of the AGK cutting rules. These rules relate combinatorial factors associated 

with different cuts of an amplitude. Applying AGK cutting rules to the system 

of four reggeized gluons we obtain the combinatorial factor of single multiplicity 

diagram(there are two such diagrams). The factor associated with every of them is 
__ j_ 

3! ' 

Summing up all contributions, using relation of discontinuity in direct channel 
to the total amplitude, we obtain : 
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Figure 6.2: Considered cuts. 



A ( ^ qq (s, t = 0) = -2sA^P< a '^P a '^ b '^ - In 2 - 

4. A. Sq 

f d 2 ^ d 2 k 2 d 2 k 3 d 2 r d% K^T^^^ 

J (2tt)3 (2tt)3 (2tt)3 (2vr)3 (2tt) 3 k?k 2 (k 2 -r) 2 (k' 2 + r) 2 ^ 

0a"xa" 2 (A.4) 

The combinatorial factor of system where two pairs of gluons form bound states 
is easy to find. We have three possibilities of pairing two gluons to form bound 
states out of four gluons. This yields the factor 2/2!. In this configuration we 
have a pomeron loop topology. In order to indicate the reflection symmetry of 
the loop we decide to keep the combinatorial factor 1/2! and to redefine vertex 
by distributing the factor 2 and also the factor of it 2 . Adding more rungs and 
including virtual corrections results in generating more logarithms of energy which 
after exponentiation can be replaced by BFKL Green's functions. Inclusion of other 
contributing diagrams that mediate transition from two glon state to four gluons 
results in replacing the kernel 2^4 with Triple Pomeron Vertex [22]. The TPV in 
our case differs slightly from the one in [22] and also in [21]. The difference is due 
to because additional color indices for two gluons coming into the vertex (when we 
consider merging) and also we included factor of n and \pl in a normalization of the 
TPV. Our motivation was to treat the vertex as locally as it was possible and that 
caused the redefinition. 
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B Angular averaging of the Triple Pomeron Ver- 
tex. 



In the forward direction the TPV takes simple form and is expressed in terms of four 
G functions. After averaging over angles, each of them gives the same contribution, 
so, to average the vertex over angles it is enough to do that for one G function. Let 
us begin with the real part of the G function 

k 2 l 2 k 2 m 2 (1 + m) 2 

Gl(1 ' m) " (k~^ + ( k + m)2 " (k-l) 2 (k + b) 2 (R1) 

Let us denote the first term in formula above A, the second B and the third C. 
An angle between 1 and m is a, an angle between m and k is (3 and an angle between 
1 and k is a+3. 
Using: 

2n 1 2tt 



o 



da — I = i 2 m ( B - 2 ) 

a + ocosa \a l — er 



we obtain: 



and similarly for B: 



,2tt f 2n ^Vl 2 

I A = I da I d(3A = — 2 T (B.3) 



o Jo 



r 



r 2w r 2w 4ir 2 k 2 m 2 
I B = da d/?g = , *™, (B.4) 
Jo Jo |m 2 - k | 

To compute integral in C = C\ + C 2 + C 3 we split it into three integrals. 



l c 2 



2w i2 



r 

dadf3—^2 : : 2 2 — (B-5) 



o 



(k 2 


— 21kcos(a 


+ (3) + l 2 ) (k 2 + 2mk cos (3 + m 2 ) 






l 2 m 2 cos a 


(k 2 


— 21kcos(a 


+ 13)+ l 2 ) (k 2 + 2mk cos (3 + m 2 ) 






m 2 



Id dB lm'cosa 
~ Jo (k 2 - 21k cos(a + (3)+ l 2 ) (k 2 + 2mk cos + m 2 ) 

f 2n m 2 
c 3 - ^ dad(3— 2 _ 21kcog ^ + p) + i 2 )(k 2 + 2mkcos/3 + m 2 ) ' b ' ' ' 



Jc z = 

Let us begin with (B.6). To compute this integral we use: 



I 



~ r , cos a 2ixb , 

da : = . B.8 

a + bcosa + csma fe2 + c 2 _ a ( a + (6 _ a ) y^pg?) 



In our case a = 1, b = — prp ; c = prp sin(/3). Inserting this in (B.8) and using 
formula which is written below to compute the integral over second angle 



Jo 



* dB ^ a = 2 - ^ (B.9) 



1 + acos ^ (i- fl )(l + a+^g) 
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we obtain: 

32/ 2 m 2 k 6 

Ic ' 2 ~ (l 2 + k 2 )(m 2 + k 2 )(l 2 + k 2 + |1 2 - k 2 |)(m 2 + k 2 ) + |m 2 - k 2 | ( } 
Applying twice (for each angle) B.2 we obtain for I Cl 

47T 2 l 2 k 4 m . 

J ci = tt2 — rrri — 5 — r^r ( B - n ) 
|1 — k ||m 2 — k I 

47T 2 m 2 k 4 

J c- 3 - TT2 — rrri — 5 — r^r l^- 1 ^) 
|1 — k ||m 2 — k I 

The total contribution is given by summing up I a, Ib, Id, Ic 2 > Ic 3 The result can 
by greatly simplified if we notice that we can consider situation when for instance 
k 2 > l 2 , m 2 . Using that condition we may drop the absolute value sign. Adding all 
terms we obtain: 

^ 21 2 m 2 k 2 - 21 2 k 4 - 2b 2 k 4 + 2k 6 

£/ = (f-k'Km.-k") (B13) 

Expanding the denominator and simplifying we obtain: 

/ = 2k 2 (B.14) 

a,..,c 3 

In case when k 2 < l 2 , m 2 we have for I Cs 

21 2 m 2 k 2 t , 

J - = (k 2 -l 2 )(k 2 - m 2 ) (R15) 

which has opposite sign to the sum of other contributions so, the total sum gives 
zero. The same situation holds in the remaining cases m 2 >k 2 >l 2 , 1 2 > k 2 >m 2 . This 
result can be simply written as: 

2?r 

.2/3/12 i r 2-s Q f m 2 1 2\ 



dadfE^l, m) = 21? 6 (V - k z )9(m 2 - k 2 ), (B.16) 



where the factor l/(27r) 2 comes from averaging. Let us now perform angular aver- 
aging of disconnected pieces of the G(l, m) function. 

I 2 m 2 

G 2 (l, m) = -k 4 In — -5(l 2 - k 2 ) - k 4 In — -5(m 2 - k 2 ) (B.17) 

(l + m) 2 (l + m) 2 

To compute the integral over angles we split the logarithm: 

I D = dad{3[ln l 2 - ln(l 2 + m 2 + 21m cos a)] (B. 18) 

Jo 
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If we assume that |1| > |m| we obtain: 

I D = 47r 2 (lnl 2 -ml 2 ) =0 (B.19) 
In case when |m|>|l| we obtain: 

l 2 

I D = 47r 2 (ln l 2 - In m 2 ) = In — (B.20) 

rri- 

Introducing 6{z) function we can write this result in a more compact form: 

I D = 4tt 2 ln(l 2 /m 2 )#(m 2 - l 2 ) (B.21) 

And similarly in the second term: 

I E = 4tt 2 ln(m 2 /l 2 )^(l 2 - m 2 ) (B.22) 

The integrated over angles and averaged G2 function reads: 

l 2 m 2 

G 2 (\, m) = -k 4 In —;d(™ 2 ~ 1 2 )^(1 2 - k 2 ) - k 4 In ^^(l 2 - m 2 )5(m 2 - k 2 ) (B.23) 
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